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ABSTRACT 

The  problem  of  scattering  of  a  scalar  plane  wave  by  a  spheroid  of  revolution 
is  solved  for  either  Dirichlet  or  Neumann  boundary  conditions,  arbitrary  major  to 
minor  axis  ratio,  and  arbitrary  incident  direction.  The  solution  is  obtained  by  using 
an  iterative  method  applied  to  solutions  of  the  corresponding  potential  problem  and 
is  expressed  as  a  series  of  products  of  Legendre  and  trigonometric  functions,  and 
ascending  powers  of  wave  number.  A  recursion  relation  for  the  coefficients  in  this 
series  is  derived.  These  results  are  employed  to  calculate  the  scattering  cross  sec¬ 
tions  for  2:1,  5:1,  and  10:1  prolate  spheroids. 
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I 

INTRODUCTION 

This  report  presents  the  complete  low  frequency  expansion  of  the  field  scat¬ 
tered  when  a  scalar  plane  wave  is  incident  from  an  arbitrary  direction  on  a  spheroid 
of  revolution  (prolate,  oblate,  or  disc),  on  which  either  Dirichlet  or  Neumann  boun¬ 
dary  conditions  are  imposed.  The  expressions  for  the  field  are  valid  everywhere  in 
space  and  for  all  values  of  the  ratio  of  spheroid  dimension  to  wavelength  within  the 
radius  of  convergence  of  the  low  frequency  expansion. 

The  work  began  as  a  demonstration  of  the  efficacy  of  a  recently  derived  tech¬ 
nique  for  solving  boundary  value  problems  for  the  Helmholtz  equation  by  iterating 
the  Green's  function  for  Laplace's  equation.  This  new  method  had  been  applied  to 
the  problem  of  scattering  by  a  sphere  both  for  a  Dirichlet  boundary  condition  (Klein- 
man,  1965)  and  a  Neumann  boundary  condition  (Ar  and  Kleinman,  1966).  The  prolate 
spheroid  was  selected  to  provide  a  more  substantial  test  of  these  methods,  which 
proved  to  work  even  better  than  anticipated. 

The  problem  of  scalar  scattering  by  a  prolate  spheroid  for  both  Dirichlet  and 
Neumann  boundary  conditions  has  been  extensively  treated.  F.B.  Sleator  (1964)  pre¬ 
sents  an  exhaustive  bibliography.  Exact  solutions  are  known  in  terms  of  spheroidal 
wave  functions  and  both  low  and  high  frequency  approximations  have  been  found.  The 
standard  methods  for  obtaining  low  frequency  approximations,  either  by  direct  ex¬ 
pansion  of  the  terms  of  the  spheroidal  function  series  in  powers  of  wave  number  or 
by  determining  each  term  in  the  expansion  as  the  solution  of  a  potential  problem 
(cf.  Noble,  i%2),  are  somewhat  cumbersome.  One  may  question  the  purpose  of 
finding  low  frequency  expansions  if  the  exact  solution  is  known.  The  answer  lies  in 
the  complexity  of  the  spheroidal  functions  which  make  analysis  and  computation  dif¬ 
ficult. 

The  present  approach,  although  certainly  not  a  trivial  calculation,  avoids 
entirely  the  use  of  spheroidal  functions  on  the  one  hand  and,  on  the  other,  obviates 
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the  need  for  solving  more  than  one  potential  problem.  The  solution  is  found  in  the 
form  of  a  series  of  products  of  spheroidal  potential  functions,  i.e.  Legendre  func¬ 
tions,  whose  coefficients  are  determinated  iteratively.  While  this  ir.  itself  might  be 
ample  justification  for  presenting  the  results,  their  value  is  considerably  enhanced 
by  the  fact  that  a  recurrence  relation  for  the  coefficients  is  found.  This  means,  in 
effect,  that  the  iteration  process  may  be  carried  out  completely  and  the  complete  low 
frequency  expansion  obtained. 

This  is  carried  out  explicitly  for  a  plane  wave  incident  from  an  arbitrary 
direction  on  a  prolate  spheroid  for  both  Dirichlet  and  Neumann  boundary  conditions. 
In  addition  to  expressions  for  the  field  valid  everywhere  in  space,  the  simplifications 
occurring  in  the  limiting  cases  of  far  zone  and  nose-on  incidence  are  explicitly  given 
as  is  the  expression  for  scattering  cross  section.  The  corresponding  results  for  an 
oblate  spheroid  and  the  important  limiting  case,  the  disc,  may  be  obtained  by  a  sim¬ 
ple  transformation  and  these  results  are  also  presented  explicitly.  Some  numerical 
calculations  of  scattering  cross  sections  of  prolate  spheroids  have  been  carried  out. 
These  results  are  presented  and  compared,  where  possible,  with  existing  data. 

In  Section  II,  the  iteration  method  is  adapted  to  take  advantage  of  the  symme¬ 
try  of  prolate  spheroid  geometry.  The  method  is  applied  to  the  Dirichlet  problem 
for  the  prolate  spheroid  in  Section  HI  and  the  Neumann  problem  in  Section  IV.  Sec¬ 
tion  V  contains  the  detailed  analytic  results  for  oblate  spheroids  and  discs.  The 
numerical  calculations  for  prolate  spheroids  are  presented  in  Section  VI.  Much  of 
the  detailed  mathematical  analysis  has  been  relegated  to  a  series  of  appendices  in 
the  hope  of  making  the  method  and  the  results  more  accessible. 

This  work  was  supported  by  the  Air  Force  Cambridge  Research  Laboratories 
under  contract  AF  19(628)-4328  and  by  the  National  Science  Foundation  under  Grant 
No.  GP  6140. 
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II 

GENERAL  CONSIDERATIONS 

In  this  section  we  present  the  problem,  the  method  of  attack,  and  some  def¬ 
initions  essential  to  a  clear  understanding  of  the  procedures  followed. 

The  problem  we  are  concerned  with  is  the  determination  of  the  scattered 
field  which  results  when  a  plane  wave  of  arbitrary  incidence  impinges  upon  a  pro¬ 
late  spheroid.  With  respect  to  a  rectangular  system  of  coordinates  (x,  y,  z),  the 
prolate  spheroid  is  oriented  with  its  axis  of  revolution  (major  axis  2a)  coinciding 
with  the  z-axis,  and  its  geometrical  center  at  the  origin.  The  minor  axis  is  2b. 
Then  the  relations  between  prolate  spheroidal  coordinates  (£,rj,  0)  and  rectangular 
coordinates  (x,  y,  z)  are+ 


X  =  c  ^(?2-  D(l-rj^)  COS0 

(2.1) 

y  =  c /(?2- lKl -r)2)  sin0 

(2.2) 

z  =  c?rj 

(2.3) 

where  c  is  half  the  interfocal  distance  of  the  spheroid,  and  1<:  £<a>,  -1  <n<+l, 
0< 0 < 2tt .  The  surfaces  i-  =  constant  represent  confocal  prolate  spheroids.  The 
metric  coefficients  of  the  spheroidal  variables  are  given  by  by 

;  h^  =  c  /(52-l)(l-rj2)  .  (2.4) 

Having  defined  the  prolate  spheroid,  we  now  turn  to  the  definition  of  the  incident 
plane  wave.  Without  loss  of  generality,  we  take  the  x-z  plane  as  the  plane  of  inci¬ 
dence.  The  direction  of  propagation  forms  an  angle  with  the  positive  z-axis, 

+  For  a  detailed  discussion  of  the  geometry  of  the  prolate  spheroid  see  Sleator  (1964). 


h5=c 


.2  2 
S--JL 

?2-l 


h  =  c 

n 


i-,2 
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(see  Fig.  1-1),  and,  if  p  is  the  observation  point  with  coordinates  (r,  0,  0),  we 
write 


i,  .  -ikrcosO 
u  (p)  =  e 


where  u  (p)  denotes  the  plane  wave  as  observed  at  p  and 


cos©  =  cos0cos0  +sin0sin0  cos0 
o  o 


(2.5) 


(2.6) 


When  0q  is  reduced  to  zero,  the  plane  wave  is  seen  to  propagate  along  the  negative 
z-axis.  The  time  dependence  is  harmonic  (e  1U^). 

We  now  state  the  problem: 

Let  S  designate  the  surface  of  a  prolate  spheroid  with  surface  coordinate 

?  ,  and  let  V  be  the  volume  exterior  to  it.  Designate  by  V  the  union  of  S  and  V: 
s 

_ _  g 

V  =  SUV.  Finally,  let  u  (p)  be  the  resulting  scattered  field  due  to  the  presence  of 
the  prolate  spheroid.  We  wish  to  determine  a  function  u(p)  such  that 
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(i) 

u(p)  =  u1(pi  +  uS(p)  , 

pe  V 

(ii) 

(V2  +  k2)uS(p)  =  0  , 

pe  V 

(iii) 

llm  r(tr'ikuS) 

r  — *■  oo  v  ' 

=  0 

(iv) 

Either 

(a)  u(pg)  =  0  , 

p  eS 
s 

or 

9u(p  ) 

(2.7) 

(2.8) 


=  0  , 


p  eS 
s 


(2. lOal 


(2. 10b 1 


Equation  (2.9)  implies  a  suppressed  time  harmonic  dependence  e  .  Moreover, 
boundary  condition  (2.10a)  refers  to  the  Dirichlet  problem  and  (2.10b)  to  the  Neu¬ 
mann  problem  and  the  two  problems  will  be  treated  separately. 

The  approach  employed  in  solving  the  problem  is  based  on  a  new  method  of 
finding  iterative  solutions  of  the  Helmholtz  equation  (Kleinman,  1965;  Ar  and  Klein- 
man,  1966).  Inherent  to  this  method  is  the  assumption  of  long  wavelength  compared 
to  the  dimensions  of  the  scatterer.  The  original  iteration  scheme  was  phrased  in 
spherical  coordinates  and  much  of  the  analysis  depended  upon  expansions  in  these 
variables.  Here  we  essentially  rederive  these  results  in  prolate  spheroidal  coor¬ 
dinates  in  which  form  the  iteration  becomes  more  tractable. 

We  start  with  a  representation  theorem  (Kleinman,  1965;  Ar  and  Kleinman, 

1966): 

Theorem:  Any  function  u(p),  defined  for  all  pe  V,  which  is  twice  differentiable  in 
V,  and  regular  at  infinity  satisfies  the  integral  equation 
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utPj)  = 


Go<pr  pv 


pv)dvt\“<ps)iGo(prps)ds 


(2.11a) 


where  is  the  normalized  static  Green's  function  of  the  first  kind,  and  the  inte¬ 
gral  equation 


u(pi)  =  \  G0N<prpv)^pv>dV 


Go<prps)to“<ps)ds  • 


(2.11b) 


where  G  is  the  normalized  static  Green's  function  of  the  second  kind.  The  normal 
o 

ized  static  Green's  function  GQ(pj,  p)  of  either  kind  is  defined  as  follows: 


(i)  V^^.p)  =  Mpjp)  , 


p^peV 


(2.12) 


(ii)  G^tp^.p)  regular  at  infinity 

(iii)  (a)  G°(p  p  )  =  0  (first  kind) 

o  1  s 

(b)  —  GN(p  ,  p  )  =  0  (second  kind)  . 
on  o  I  s 

The  normal  is  directed  out  of  the  volume  V.  Moreover,  we  define  a  function  f(p)  to 
be  regular  at  infinity  if  it  satisfies  the  Kellog  (1929)  conditions 


lim  |rf(p)j  <  oo  and 
r  — >oo 


lim 
r  — >oo 


r 


2  §fl£) 
9r 


<  oo  , 


0 

2n  . 


(2.13) 


Using  expressions  (2.4),  it  can  be  readily  shown  that  in  prolate  spheroidal  coor¬ 
dinates 
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poo  P+1  r\2it 
uXPl)  = -2lkc2  \  A  drA  d0G°(pi,p)^2-l)^-); 

J?  J-i  Jo 


2 

or) 


>+l  p 2^ 


dr)  ^  dMPs)  G°(pr  ps). 


+  n)u(p)J 


(2.20a) 


for  the  Dirichlet  case,  and 


^Pj)  = 


>oo  p+1  p2ff 

-2ikc2\  d?  j  dr)^  d0  G^(Pl,p)|^ 


..2  aUp)  _  2  3u(p)  , 

(?  +  (r)  -  1)  +  (5  +  n)u(p) 


'+1  P27T 

♦<»\  M  d*Go(prPs)l|'“(ps) 


(2.20b) 


'-1  Jo 


s 


for  the  Neumann  case. 

These  are  the  integrodifferential  equations  that  we  have  to  solve.  The  first 
one  involves  the  normalized  static  Green's  function  of  the  first  kind  (Dirichlet  boun¬ 
dary  condition)  defined  by  (2. 12)  and  given  by 

Gr<Pl-p)=-lh£X](-1)m£m(2n  +  1)[(tTS)0 

n=0  m=0  L- 


T,m.  .  in,  . 
P  (n,)P  (n) 
n  1  n 


*~r  Pm(€,)QmCC)'  , 

n  i  n  |  pm(5  } 

Qm(?.)Qm(?) 
m(„  t  n  1  n 


l^p 

n  n  1 


Q“  (?  ) 

n  s 


€>€, 


€<€, 


(2.21) 


|iee  for  example,  Morse  and  Feshbach  (1953,  p.  1291).  The  existing  differences  are 
due  to  a  different  normalization  and  a  different  definition  of  the  Legendre  functions. 


8 


THE  UNIVERSITY  OF  MICHIGAN 

7133-5-T 


The  corresponding  Green's  function  of  the  second  kind  is  of  similar  form  except  for 
involving  the  ratio  of  the  derivatives  of  Legendre  functions  so  that  the  boundary  con¬ 
dition  (2. 12. iiib)  is  satisfied^  Equation  (2.20b)  involves  the  normalized  static 
Green's  function  of  the  second  kind  (Neumann  boundary  condition)  defined  by  (2.  12) 
and  given  by 


ao  n 


GX-»,= 


T'~~ ’(-l)me  <2n+l)  cosm(0  -0) 


„m,  ,_m,  . 

p  (nJPtn) 

n  1  n 


- 


*—  n  n  l 


Pm(C  )' 

Q  (?  )'  n  1  ° 

n  s 


€  >£ 


l 


,  (2.22) 


where  a  prime  on  a  function  denotes  differentiation  with  respect  to  The  symbol 

€  is  the  Neumann  factor  defined  by 
m 


c  =  A 
m 


1 ,  m  =  0 

L  2,  m  =  1,  2,  3, . . . 

The  associated  Legendre  functions  are  defined  as  follows: 


(2.23) 


pin.  .  _  _1 _ Hn+m  +  1) 

n  ^  2m  fTn- m  +  l)f(m+  1) 


(l-^2)m^2  2F1(m-n,  n  +  m  + 1;  m+1;  ) , 


m  n  r<n+i> 

PnM-2  n^T 


iw-ll  <2  .  (2.24) 

27  ,  2  .vm/2  n-m  _  f  m  -  n+ 1  m-n  1  1  \ 

viraji)  w  _1)  “  2FiV“~T—  • 

|m|>1;  t  arg<^±  1 )  j  <  tt  .  (2.25) 


i 
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m.  _  (-l)m  f(n+m-H)r(l/2)  (u-  l)m^2  ( o+m+2  n+m+1  3  A 

n  _n+l  r.  3..  n+m+1  2  IV  2  '  2  ;n  2: 

2  l(n+-)  M  M 

j#i|>l;  |arg<^ -1)|  <  .  (2.26) 


Definitions  (2.25)  and  (2. 26)  agree  with  those  given  by  Magnus  and  Oberhettinger 
(1949,  pp  64  and  60,  respectively),  while  (2.24)  differs  by  a  factor  of  (-l)m. 

To  solve  the  integrodifferential  equations  (2.20)  for  Up)  we  proceed  as  fol¬ 
lows: 

We  write  Up)  as  a  power  series  in  k  of  the  form 


ikc)Mu  (p) 
M 


and  we  substitute  in  equations  (2. 20)  to  obtain  an  iteration  scheme  for  u>  (p).  We 

M 

subsequently  show  that  these  coefficients  of  k  are  of  a  particular  form  and  develop 

recurrence  relations  through  which  u  (p)  can  be  found  for  arbitrary  M. 

M 
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THE  DIRICHLET  PROBLEM 

3.1  The  Iteration  Scheme 

The  appropriate  integrodifferential  equation  for  the  Dirichlet  problem  is 
(2. 20a)  which  we  repeat  here  for  convenience 


u(p^)  =  -2ikc 


cd  r+1  r>2  it 


-1 


d?\  drj \  djDG^Pj.p)  (?2-  1) +  (q2- 1)^^ qMp) 


+1  rh 


-<»)  dr^  dMp8)^°o<p:'p8)- 


(3.1) 


The  appropriate  Green's  function  is  given  by  (2.21)  and  the  boundary  condition  satis¬ 
fied  by  u(p)  is  seen  to  be,  from  equations  (2.7),  (2.10a)  and  (2.14), 


-ikc(?  +q) 

u(ps)  =  -u  (pg)e 


(3.2) 


The  incident  plane  wave  u  (p)  is  given  by  (2.5)  which  can  be  written  in  prolate  spher¬ 


oidal  coordinates  as 


i.  v  -ikrcos0 
u  (p)  =  e  =  e 


-ike  [cos  %r) +  sin  Ji  - n2  cos0] 


(3.3) 


Denote  the  surface  integral  of  (3. 1)  by  I  (p^): 


p  +  1  p27T 

v  __  _,.2  v\  \  j  .  %  d  „D. 


-1  JO 


KPl)  =  -c(?s-i)\  dn\  d0Upg)  — Go(Pl,ps)  . 


(3.4) 


In  Appendix  B  we  show  that  I  (p^)  may  be  written  in  the  following  form 
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-ike?  v-"1 

IS<P1)=e  8  /  ,  (-ike)  ^(Pj) 


(3.5) 


where 


M  i 


^Pl*  =  Z_.Z_o  )p{%  )Q™(5  )cosm0 

1=0  m=0  811,1  1 


(3.6) 


A“,m(5  )  - 

l  S 


r(?s-c°s9o>  rr 
em^  0M+1  +  (f+m)'. 


A  cos  9  -  1 
,m  /  s  o 


f  V  ?  i  COS0 


M+f  even 


0  ,  M  +  i  odd  . 


(3.7) 


Note  that  1^  is  independent  of  k.  Moreover,  let 


+ick? 

0(  Pj )  =  e  su(p  ) 


(3.8) 


where  0(Pj)  is  assumed  to  have  a  power  series  expansion  in  k  of  the  form 


03 

■  £<-. 


ikc)M^M(p)  . 


(3.9) 


Substitution  of  (3.9)  in  (3.8)  and  the  resulting  equation  together  with  (3.5)  in  (3. 1) 


gives 
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The  interchange  in  differentiation  and  summation,  and  summation  and  integration 
was  made  by  assuming  (3.9)  to  converge  absolutely  and  uniformly  and  to  be  term  by 
term  differentiable  with  respect  to  each  of  the  variables  and  the  resulting  series  to 
be  uniformly  and  absolutely  convergent.  Collecting  the  coefficients  of  equal  powers 
of  k  in  the  above  equation,  we  arrive  at  the  following  iteration  scheme: 


w = 0pi> 


(3.10a) 


■\00  r>+l  r\2jr 


,D, 


tf'M+l^l*  =  2c  \  d?  \  dr5  \  d0  G0<Pr  P> 

■1  Jo 


+(? 


(C2  .)aVp>  2  3Vp) 

‘5  +  (p 

+  1M+1<P1>  •  (3.10b) 


3.2  The  Recurrence  Relations 

We  shall  now  solve  for  the  M+  1st  iterate  in  (3. 10b).  In  order  to  do  this  we 
need  to  establish  the  fact  that  &  may  be  written  as 


M  M  t 


Vp)  = 


t=0  r=0  1=0 


D^'/(?  )Qf(C)pJ(rj)cosl0 
r,t  s  r  i 


(3.11) 


for  all  M  (M  =  0, 1,  2, . . . ).  This  is  accomplished  using  the  principle  of  mathematical 
induction,  that  is,  first  we  show  that  (3.11)  holds  for  M  -  0  and  secondly  we  show 
that  if  it  holds  for  M,  it  also  holds  for  M+l. 
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That  the  representation  holds  for  M  =  0  is  obvious  since,  with  equations 
(3.10a)  and  (3.6), 


<M P)  =  A°*°(?  )Q  (0 
o  o  s  o 


(3.12) 


which  clearly  is  of  the  form  (3.11). 


Next  assume  that  (3.11)  holds  for  M.  We  wish  to  show  that  may 


then  be  written  as 


M+l  M+l  t 


W1* 


-  )Q'<?)p'(n)cosif> 

i_A  _ A  |_A  *  »  1  OX  l 


t=0  r=0  1=0 


The  analysis  which  establishes  this  is  somewhat  tedious;  however,  in  the  process 

we  actually  arrive  at  an  expression  for  D®**1,1  in  terms  of  DM,i  which  in  fact  is 

r,  t  r, t 

the  major  goal  of  this  section. 

First  note  that  the  second  term  in  (3.10b)  has  already  been  shown  to  be  of  the 
form  (3. 11)  jsee  (3.6),  (3.7)J  .  Next  denote  the  volume  integral  of  (3. 10b)  by 

v 

I^+1(p,)  and  substitute  in  it  the  Green's  function  of  (2.21).  Then, 


aa  n 


M+l^l*  "  2n  L— j( 

n=0  m=0 


•Dmc  (2n+i)Ra^n2pm( „,) 
m  L(n+m):J  n  T 


pOO  p+1  p2ff  _ 

\  dCC^(C,Ci,Cs)  \  drjP™(n)\  d^cosm(0-01>  (?2-  D 

J?  J-l  Jo 


2  ^M(P) 

a? 


2  a/M(p) 

+  '1)~^r_  +  (?  +  n)^M(p) 


(3.13) 
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)  = 

n  Is 


fp”(?l)Q”<5) 


J 


Pm(?  )  > 

zr  ”  q!X>Q“<?> 

<C(V  *<« 


Substitution  of  (3. 11)  in  (3. 13)  leads  to 


co  n  M  M  t 


f~(n-m)rf 

L(n+m)!  j 


IM+1(P1>  "  2jt  Z_i  Z_iZ_.  Z_iZ^,(_1)  em(2nfl)[ 

n=0  m=0  t=0  r=0  i=o  L 

•DM,5s,P>1>^  d5C<5’5 l-V^  d”P>> 


"  o  -  dQ^(0  o  i  dP 

d0 cos  m(0  -  0  )  cos  f  0  (?  -  l)P  (n) . -  +(rj  -1)Q  (?)  — ~ 

i  t  d?  r  d 


+  (?  +  n)Qr(?)P^(n) 


Performing  the  angular  integration  and  rearranging  terms  we  get 


oo  M  M  t 


*M+1*P1* 


n=0  t=0  r=0 


poo  p-fi 

pVjcosfjJ,  \  d?c'(?,f,U  dnP^n)f|V-l)4+fl 


THE  UNIVERSITY  OF  MICHIGAN 

7133-5-T 


To  perform  the  r)  integration  we  use  the  relation  (Magnus  and  Oberhettinger,  1949, 


pp  61-62) 


(n+lXn-l  +  1) 
2n+ 1 


n+1 


^Qn+1 


n(n+l) 
2n+  1 


rV  (z) 
n-1 


A-i(z) 


n,l  =  0,1,2, 


(3.16) 


Substituting  this  relation  in  (3.15),  we  get 
oo  M  M  t 


oo 


«+l 


d5  \  dn  p'<h)  \  4iif^  <,<n><(5> 


1-1 


oo  M+l  M  t 


=  (2n+1>  fefT)  P^n^cosf^ 

n=0  t=0  r=0  i=0  LV  ,,J 


oo 


<+l 


J, 


d?C  (?,?.,?  ) 
n  1  s 


dTjp'(n)pJ(n)  \  h 

■  l  *- 


2t-  1  r,  t-r  s  r 


+  DrV({U 

r,  t  s 


«2-  +?J  q'<?>  t  li±^r±i2  d“;11<cs)q'<5>! 


(3.17) 
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M  f 

where  in  this  last  expression  we  have  adopted  the  convention  that  D  '  (?)  is  iden- 

r,  t  s 

tically  zero  whenever  any  of  the  subscripts  or  the  other  superscript  is  greater  than 
M.  We  now  employ  the  following  orthogonality  property  for  the  Legendre  functions 
of  the  first  kind  (Magnus  and  Oberhettinger,  1949,  p.  54), 


'+1 


.  J,  ,  2  (n+f)i  . 

dxP  (x)P  (x)  =  — rr  6 

n  m  2n+l  (n-/)i  nm 


1-1 


to  obtain 


M+l  M  t 

IM+i(pi)  =  ~2  E  E  XE  (t+jj.  pt<n1)cosi01 

t=0  r=0  1=0 


dS  c'(f  f  5  )  fl T-  dM'  '  (f  )  +  dM.  I  .  T1  «(5) 

®  2t  - 1  r.t-l'V'  2t+3  r,  t+l''s^r 

Employing  once  more  the  relation  (3.16)  in  the  equation  above,  we  write 
M+l  M  t  poo 

‘m+i(>V  ■  -2  5  S  oTTTt  \  V 

t=u  r=0  i=u 

s 


-  *<t -J1  I:  l 

+  ITl  Dr,t-l(V 


-  2t  +  3  Ur,t+l(V  Qr  *  2r+l  Dr,t  (VQr-l(5)1 


(cont'd) 
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M+l  M+l  t 


•  A  ■  A  V 

V  VY"1,  ni  (t-l)'  I,  , 

I -f  u^n  ?t(VC08,M  l^rr  Vi.  t<5f 

t=0  r=0  1=0  L 


t(t-D  M,l  +  (t+lHt  +  l  +  U  nM,l 
+  2t-l  r.t-rs*  2t+3  r.tfl'V 


•lMr.+l  +  l)  M.l  ) 
2r  +  3  i+l.t'Vj 


i  i  s  r 


(3.19) 


In  arriving  at  this  last  expression,  one  must  bear  in  mind  that  ^  (?g)  is  iden¬ 
tically  zero  whenever  r,  t  or  <  is  greater  than  M.  As  shown  in  Appendix  C, 


oo  p  ^  - 

-rtr,i)V„  Hi  • 

€  LW 

s 


r  i  t 


(3.20) 


Furthermore,  whenever  r  =  t  in  (3.19),  the  bracketed  coefficient  is  equal  to  zero. 


This  follows  from  the  fact,  established  in  Appendix  D,  that  the  relationship 

DM,i(?  )  =  (tl)rHDM’<(?  ) 
r,  t  s  t,  r  s 


(3.21) 


holds  among  these  coefficients.  Thus  we  need  not  evaluate  terms  in  (3.19)  when 
r  =  t.  Substituting,  then,  (3.20)  in  (3. 19)  we  get 
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M+l  M+l,  t 

IM+l<pl)  =  *2  S  ^  rtr+D- 
t=0  r=0  x=0 


t(t+l) 


2r- 1  r-l.t'V 


_  M,  I  ,  ,  (t+lMt  +  i  +  1)  _M.i 

2t- 1  r,t- 

(r+lXr+f +1)  M,  I 


_  ra,  *  ,  _  .  +  -  .  . 

+  2t-l  r,t-l  V"  2t  +  3  r.t+l'V 


2r+3 


Dm.,‘V 


i._  .  W'W 


L<V5s> 


Pt(n1)cosi01 .. 


(3.22) 

where  the  prime  on  the  summation  for  r  indicates  that  the  term  t  =  r  must  be  de¬ 
leted.  This  may  be  rewritten  as 


M+l  M+l  t 


M+l.  I, 


w»i>  ■  4 44FrV  ■ 


(3.23) 


whe.-e 


M+l,f  2  [~ rir-t)  nM,f  _  t(t-i) 

r,t  s  r(r+l)-t(t  +  l)  j_2r-l  r-l,i  s*  +  2t  —  1 

+  (t  +  l)(t+i  +  1)  M.i  (r+lKr+l+1)  M.i  7] 

2t  +  3  r,  t+i  s*  ”  >r  +  3  r+l,t  s_| 


r,  t-1  s 


r  ^  t 


(3.24a! 


M+l-  Q*{E  ) 

em+1’'(5s)  -  -  J2  -p"  • 

r=«  q;<€  )  ’ 

t  s 


t.t 


(3.24b! 


and  the  prime  on  the  summation  indicates  that  the  term  r  =  t  must  be  deleted  (see 
Appendix  E).  Thus  (3.10b)  can  be  written  as 
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M+l  M+l  t 


WP1>:=  WP1)+WP1 


)  = 


t=0  r=0  f=0 


V(?)pJ(n)cos^ 

i|  i  s  r  i 


(3.25) 


where 


dm«.«  =  eM|U 
r, .  r,  t  s 


=  em+1,((5  )+am+1'' 

t,t  S  t 


r  =  t 


Equation  .25)  is  clearly  of  the  form  (3.11)  which  is  what  we  wished  to  establish. 

Not  only  have  we  completed  this  inductive  proof  but,  in  the  process,  we  have  derived 

M  I 

recurrence  relations  for  the  coefficients  D  ’  (£  ): 

r,  t  s 


\ _  f ijr-l)  M ,/  (p  .  _  t(t-OnM,f 

-t(t+ 1)  |_2r-  1  r-1,  is  +  2t  —  1  r.t-l'V 


_  M+i ,  f . .  ,  _ ; 

Dr,t  r(r+ 1) 


+  (t  +  l)(t  +  <  + 1)  M,<  (r+l)(r+f  +  1)  M.f 

2t  +  3  r,t+l'V"  2r+  3  r+l,t'V 


r  i  t 

M  =  0,1,2,  ... 


(3.26a) 


DM+1-'(f)=-E 


t,t 


M+1  •  Q^(?  ) 

W  „M+l,f . M+l,  f  %  ..  „  .  „ 

— - D  '  (?  )  +  A  ’(£);  M  =  0,1,2,... 

i .  _  .  r,  t  s  t  s 


r=0  Qt<y 


(3.2Gb) 


with 
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M  i 

where,  in  both  of  the  above  equations,  D  '  (?)  is  given  by  equations  (3.26a,  b,  c) 

M  t  Tlt  8 

and,  in  turn,  A  '  (?  )  is  given  by  (3.7). 

t  S 

3.3  The  Far  Field  and  the  Scattering  Cross  Section 

From  the  definition  of  Q^\?)  *n  (2.26)  and  (3,28)  the  far  field  is  given  by 


ike?  ao 


n-M  M  t 


si.  .  e 
u  (p,)= 


w  U  _  .U  ill  XV*  l 

£(-lkc)°  X]  ~  XX'-1^ 

?1  n=0  M=0  1  t=0  1=0  0,1  8 


Pt(n1)cosfjD  .  (3.29) 


Since  the  incident  wave  is  of  unit  amplitude  and  r  ~  c?  in  the  far  field,  the  scatter¬ 
ing  cross  section  is  given  by 


<j=  lim  4xrc2?2  uSf(p  )|2  =  4;r 

p-  A  1 


?^— >oo 


0D 


n=0 


(-ikc)nuSf(p  ) 
n  l 


(3.30) 


where 


n  _  .n-M  M  t 

(S  +0,) 


%t(p,>  -  X]-' IXm)'  d”'1!;  )P%  >CO3<0  . 

n  1  M=0  m  t=0  1=0  o,t  s  .  1  1 


(3.31) 


Assuming  k  real,  we  can  rewrite  (3.30)  as  follows: 
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oo 

£- 

n=0 

CD 

£<- 

n=0 

ao 

£« 

n=0 

ao 

£<- 

n=0  m=0 


2 

)°uSf(p  )  =  47rc2  />  !(-ikc)P 

n  1  /~3 

n=0 

ao 

uSf(p  )  (-ikc)°  uSf(p, ) 

n  1  ITO  n  1 

oo 

^^(p.)  y~",(-ikc)Pu^f(p) 

n  1  C - rr1  n  1 

n=0 

ao 

n  uSf^Pi  ^  y  ,ukc)p  ubf(p ) 

n  1  n  1 

n=0 

n 

l”  £  (-l)muSf  (p,)USf(n) 
c — £  n-m  1  m  1 

oo  2n 

£^£ 


=  47rc"X  ,(kc) 

n=0 


m=0 


(3.32) 


3.4  Nose- on  Incidence 

In  the  case  of  nose-on  incidence  (0  =  0)  quite  a  few  simplifications  occur. 

o 

If  we  set  0  =0  in  (3.7),  it  becomes  obvious  from  the  definition  of  the 
° 

Legendre  function  Pm(p),  I yu -  1 1  <2,  equation  (2. 24),  that  A.  ,m(?  )  becomes  zero 
n  1  1  is 

unless  m  =0.  We  then  conclude  that  in  the  case  in  which  the  incident  plane  wave 
propagates  along  the  z-axis  there  is  no  dependence  on  the  azimuthal  angle  0.  This 
simplifies  the  results  as  follows: 

Equation  (3.6)  can  be  written 


M 

■>i> = £  Ar<vpl(vQt<51)  • 


with 


(3.33) 
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( ?  -  I )  ( 2t  +  1 ) 

o 


am(?  )  = 

t  s 


^F^W)Tr(mi+i)q, 

^  0  .  M  + 1  odd 


<VV 


Equations  (3.26a,  b,  c)  become 


Dr,t(V  r(r+  1)  -  t(t  + 1)  _  2r-  1  °r-l,  t(' V  +  2^1  Dr|t-l(V 


M  +  t  even 


+  liULL  dm  <f  )  (r+l)2  M  ,  1  r*t 
2t  +  3  r,  t+1  s  2r  +  3  Dr+l,t(V  ;  M  =  i 


M  =  0,1,2,  . 


The  scattered  field,  us(Pi).  becomes 


U  (Pj )  =  e 


and  the  far  field. 


/c  -  »n-M  M  M 

W  U  fa-W!  ^4.°,t 


(WVW 


sf  e 

u  <P,>  =  — 
1  $ 


ike?  co 

1  v — I 


i  — i  i  [p  \n-M  .  M 

i  ^  ('IkC)  S  DM(Vpt(V  • 
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IV 

THE  NEUMANN  PROELEM 

4. 1  The  Iteration  Scheme 

The  appropriate  integrodifferential  equation  for  the  Neumann  problem  is 
(2.20b)  which  we  repeat  here. 


/(.2  3u(p)  _  ,  2  iV  3u(p) 

(5  + 


p+1  p 2n 

+  (f  +  nMp)j+  c(?g- i)  \  dn  \  d0G^(p1,ps)^|-u(ps)  , 

J-i  Jo  s 


f 

>oo  r 

>  +  !  (• 

>2jr 

2  1 

u(p  )  =  -2ikc 

d?  ' 

dn 

d0G*(p  p> 

V 

k  - 

)-l  V. 

Jo 

(4.1) 


,N. 


with  G^p^,  P)  given  by  (2.22).  The  appropriate  boundary  condition  is  given  by 
(2. 10b),  which  through  equations  (2.17)  and  (2. 18),  may  be  written. 


1  l5llL  JL 

c  |  .2  23? 

vn  s 


ikc(?  ±n> 

u  (p  )  +  e  S  oXp 

s  s 


n... 


(4.2) 


Excluding  the  case  in  which  the  prolate  spheroid  degenerates  to  a  wire  of  finite 
length  (?  =1),  we  can  write 


Mu‘ 


ikc(?gt  q) 

(p  )  +  e  u(p  ) 

s  s 


=  0 


(4.3) 


from  which 


3uXp  )  -ikc(?  t rj)  Ou^p  ) 

if-n1-  5  ir- 


(4.4) 


Substitution  of  (4. 4)  iri  (4.1)  leads  to 


26 


THE  UNIVERSITY  OF  MICHIGAN 

7133-5-T 


„  (?  -  cos  6  r'1 

.  M,m,  _  .  _  I — _ s _ o _ 

.  <?S  -  +  Q,X>' 


(2f+l) 


(i  -  m)! 
(f+m)’. 


/?  cos0  ti\  r 

.  MW  -S. — 2—  )  +  — L_  t 

•  V  ?  t  cos  0  )  -  2 

Ns  O  '  5  -1  L 


/?  COS©  tl 

<(?  cos 9  ±1)  pff-f-r— V- 

S  O  *  \  ?  COS0 

vs-  o 


A  cos  0  +1Y] 

-(f+m)(£  tcos©  )Pm  -~z — )  \  for  M+  f  even, 
s  o  f-l\^?slcos0  ^/j 

(4.9) 


aJ1*  m(?  )  =  0  , 

I  s 


for  M+f  odd 


(4. 1C 


From  now  on,  the  procedure  for  developing  an  iteration  scheme  parallels  that  of 
the  Dirichlet  problem.  After  writing 


^Pj)  =  e  uXp^ 


(4.11 


and  assuming  a  low  frequency  expansion  in  powers  of  k  for  iXp^), 


cXp.)  -  /  ,  (-ikc)MC  (p  )  , 

M  -  0 


(4. 12 


we  substitute  (4.11)  in  (4. 12)  and  the  resulting  expression  together  with  (4.  7)  in 
(4.5).  Equating  coefficients  of  equal  powers  of  k,  we  obtain  the  following  iteration 


scheme: 


t  (P.)  =  0 


(4.  13 
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2?r 

d*GXps)'W 

0 


M  =  0,1,2, ... 


(4.13b) 


4.2  The  Recurrence  Relations 

The  procedure  we  shall  follow  here  is  practically  identical  to  that  for  the 

Dirichlet  case.  We  assume  ^.,(p, )  to  be  of  the  form 

M  1 


M  M  t 


w  *  ZEE^VX'p!(v«^i  •  « -  o.  i .  ^ 2-  •  ■  ■ 

t=0  1=0  ' 


(4.14) 


which  we  substitute  in  (4.13b)  and  solve  for  ip...  ,(p,).  If  .(p.)  turns  out  to  be 

M+l  1  M+l  1 

of  the  form  (4. 14),  then  because  of  0  (p,)  being  zero  we  can  conclude  that  (4. 14)  is 

o  1 

true. 

The  volume  integral  of  (4. 13b)  is  practically  identical  to  that  of  (3.10b)  ex¬ 
cept  for  the  Green’s  function.  From  equations  (2.21)  and  (2.22)  we  see  that  these 

two  functions  are  identical  except  for  their  dependence  on  the  surface  coordinate  £  , 

s 

If  we  denote  the  volume  integral  of  (4. 13b)  by  Ij^+j(pj),  we  can  1186  the  result  of 
(3.19)  and  write 
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M+l  M+i  t 


‘W  *  -2  7777?  p'lt,*008'!*! 

t-0  r=0  1=0  1  ’■ 


Hr-  t)  M. f  _  t(t-<>  M, (  +  (t  +  lKt  +  l  i  1)  nM.I  .. 

2r-  1  Dr-l,t(y  ITT  Dr,t-l(V  “  Dr,t+1(V 

'cn 

(r+lKr+t  +  l)  nM.<  I  \  .  / 

2r+3  Dr*l,t(Vf  \  *  V5'?l'  cs)Qr 


(4.15) 


where 


P<(f).  «>*, 


(4.1G) 


w 


?  <c 


i 


M  I 

It  is  understood  in  (4. 15)  that  D*  '  (?  )  is  identically  zero  whenever  r,  t  or  i  is 

r,  t  s 

greater  than  M. 

From  Appendix  G, 


100 


,  *  c  c  r  \rJiz\ _ (-Df  (t  +  f); 

dcKt(^?r  VQr(=)  r(r+  1)  -  t(t+  1)  (t-f)l 


W  l 

~p-  Qt(5l)-«r(?l> 

W 


r  =f  t  .  (4.17) 

Furthermore,  through  an  inductive  argument  identical  to  that  given  in  Appendix  D 
for  the  Dirichlet  case,  w'e  can  show  that 


DM,  £(F  )  =  (+  ) 

r,  t  s  t,  r  s 


(4.18) 


Employing  (4. 17)  in  (4. 15),  we  obtain 
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M+l  M+l  .  t 


WlV  =  "2  ^-Q->  ^  r(r+ 1)  -t(t+Tj  j 


X  nM*1  (C  V  +  (t+l)(t+i  +  l)  M.l 

+  2t-l  r,t-l  V'  2t+3  r.t+1'  s 


(r+lHr+l  +  l)  M  ,* 

2r+3  rH.t'V 


V(€  )'  i  i  "1  . 

“l^*qi(€I>“Qr<€l)|Pt(nl)c08l^l 

i-Q.(C  )’  r  J  1  1  1 


(4.19) 


Having  evaluated  the  volume  integral  of  equation  (4. 13b)  we  now  turn  to  the  surface 
integral  of  the  same  equation  and  we  denote  it  by  Kp^). 


p+1  p2r 

I(P,)  =  c(5^-l)\  dn\  <0<#P1-p8WM<p8)  . 

Substituting  equations  (2.22)  and  (4.14)  in  (4.20),  we  obtain 

oo  n 


(4.20) 


I(p,)  =  -“Uf-H/,  (2n+ 1) 

1  457-  S  L - FT*  m 

n=0  m=0 

M  M  s 


(n-n)1.  2 
Jn+rj)!_ 


_m,  . 

PD  ‘V 


•  y  y  ,y  |dm'/<?  >q'<?  pm<?  >-— —  Qm<;  > 

trfW™  rA  '  r  ‘  1  1  L"  ‘  Qm(f  >■  n  SJ 

n  s 

n2ir 


X  ^  dn  p“(r})Pt<n)  ^  d0cosm(0-0j)cosf0  . 

Using  equation  (B.4)  for  the  Wronskian  ard  at  the  same  time  performing  first  the 
integration  with  respect  to  0  and  then  the  integration  with  respect  to  rj  according 
to  (3. 18),  we  obtain 
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YllYYL )cosi0 

t=0  r=0  f-0  *  Q  (f  )*  1 

t  s 


(4.21) 


From  equations  (4.19)  and  (4.21),  it  is  clear  that  0  (p  )  of  (4. 13b)  is  of  the  form 

m'l  1 

given  by  (4. 14).  At  this  point,  then,  we  not  only  have  concluded  the  inductive  argu¬ 
ment  that  the  representation  (4, 14)  of  c-  (p^)  is  correct,  but  in  exactly  the  same 
fashion  as  in  the  Dirichlet  problem  we  end  up  with  the  following  recurrence  relation¬ 
ships: 

)  _ _ 2 _  [rjr-t)  M.f  -  tft^O  M.f 

Dr,t  (V  rf  r+  1)  -  t(t  + 1)  2r- 1  Dr-l.t(V  +  2t -  1  Dr,t-l(V 

+  (t  +  l)(t  +  (  +1)  M.f  (r  +  l)(r  +  f +1)  M.f 

2t  +  3  r,t+l  s  ”  2r  +  3  rH.t  sj  ' 


M  =  0,1,2,... 

r  J  t 


(4.22a) 


^  Q1*?  )’ 


D,M:M(n  =  -E  ^ dm;^(5 )+AtM+i*v ), 

M  s  ^  q1(?  )'  r’1  5  Uq  (?)’  r'1  s  1 

t  s  t  s 


M  =  0,1,2,  ... 


(4.22b) 


D°’  V  )  =  0  , 

o,  0  s 


(4.22c) 


where  A  ’  (£  )  is  given  by  (4.9a,  b).  The  prime  on  the  first  summation  in  (4.22b) 
t  s 

denotes  that  the  term  r  - 1  must  be  deleted. 
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The  scattered  field  u'(p^)  for  the  Neumann  problem  is  given  by  the  same 

expressions  as  for  the  Diriehlel  problem  (equations  (3.27)  and  (.3.28))  with  the  under- 

M  f 

standing  of  course  that  the  coefficients  IV  (£  )  are  this  time  given  by  equations 
(3.22).  The  same  is  true  for  the  far  field  and  the  scattering  cross  section  expres¬ 
sions  (see  Section  3.3). 

4.3  Nose-on  Incidence 

When  0=0,  (4.9a)  becomes  zero  unless  m  =  0.  This  is  so  because  of  the 
o 

definition  of  the  Legendre  function  (/a)  for  |^-l|  <2,  equation  (2. 24).  Conse¬ 
quently,  when  the  incoming  plane  wave  propagates  along  the  z-axis,  there  is  no  de¬ 
pendence  of  the  scattered  field  on  the  azimuthal  angle  a  result  we  should  expect 
since  the  z-axis  is  the  axis  of  symmetry  of  the  prolate  spheroid.  Due  to  the  sub¬ 
stantial  amount  of  simplification,  we  redefine  our  results  for  the  Neumann  problem 
as  follows: 

Equation  (4.8)  becomes 


AI 


'mV  - 


M 


t=0 


A,  •Vp.Wi1 


M  -  1,2,3, 


(4.23) 


where 


AiX>  ■ 


-fir  M 


(2t  +  1 )(?  -1) 

s 


M-l 


M+l  /  M  -7 


-  AI±1  +  i  j  Q  (£  ). 

1  v  2  2 J  V  s 


M  + 1  even 


(4.24) 


M  +  t  odd. 


M 


The  scattered  field  u  (p^)  is  given  by  (3.89)  witn  D‘  ^(?  )  given  by 
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DMt!(?  )  = 


r,t  s  r(r+l)-t(t  +  l) 


'2  2 

r  _  t  A!  . 

2r-  1  Dr-l,t(y  +  2t-l°r,t-l(V 


-  2t  +  3  r.t+l'V  2r+  3  rH.t  V 


M  =  0,1,2, 
r  f  t 


(4.25a) 


M+l..  .  ^T"1  Qr  V  M+l  .  "S-1  M  .  M+l. 

D.,t  <5s’  *  -2-f  WT  Dr,.  (5s)  +  ^  qTT)7  DrA>+A.  «.>  • 
r=0  t  s  r=0  t  s 


M  =  0,1,2, ... 


(4.25b) 


with 


D°  (?  )  =  0  .  (4.25c) 

o,o  s 

The  expression  for  the  far  field  is  the  same  as  the  one  for  the  Dinchlet  problem 

M 

(equation  (3.37))  with  D  (£  )  as  above.  The  same  is  true  for  the  coefficient  of  the 

o,  t  s 

sf 

scattering  cross  section,  u  (p  ),  which  is  given  by  (3.38). 
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V 

THE  OBLATE  SPHEROID  AND  THE  DISC 

The  method  employed  in  the  preceding  sections  to  determine  the  field  scat¬ 
tered  b}  a  prolate  spheroid  can  be  employed  in  a  straightforward  manner  to  deter¬ 
mine  the  field  scattered  by  an  oblate  spheroid.  This  is  not  necessary,  however, 
since  we  can  transform  the  prolate  spheroid  into  an  oblate  spheroid  and  utilize  the 
results  already  obtained  to  determine  the  field  scattered  by  the  oblate  spheroid. 
Specifically,  if  we  let  £  — *  i£  and  c— >-ic,  the  prolate  spheroid  is  transformed  into 
an  oblate  one  with  the  axis  of  revolution  (minor  axis  2b)  coincident  with  the  z-axis  of 
a  rectangular  coordinate  systems  (cf.  Morse  and  Feshbach,  1953,  p.  1502).  The 
ranges  of  the  variables  now  are  0  ^  <  oo,  -K<  n  £  +1,  and  0  £  0  £  2n  .  Moreover, 


FIG.  5-1. 


if  we  let  *  0,  the  oblate  spheroid  degenerates  into  a  disc  of  infinitesimal  thick¬ 
ness,  radius  c  (the  semifocal  distance),  and  coplanar  with  the  x-y  plane.  In  the 
remainder  of  this  section  we  shall  treat  each  body  separately. 
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5. 1  The  Oblate  Spheroid 


With  the  incident  field  given  by  (see  equation  (3.  3)), 


i,  .  -Her cos  0 

u  (p)  =  e  -  e 


-ike  cos  0 
L  o 


£r;+  sin0(  >r«  -1  z1  -  T)  COS  j 


the  scattered  field  is  given  by  (3.28)  with  ?  =  i £  and  c  =  -ic  and  can  be  written  as 


follows 


_  to  n 

T,  •„  •»  V,  -v 

u  (p  )  =  e  /  (-ike;  /  ,(-i) 

f  n  «  •  A 


n  _  .n-M  M  M  t 

y\  vvv.»i.y. . 

fes1  r-'  s 


Qr(iV?t(VCOS^l  * 


(5. 1) 


where,  for  the  Dirichlet  problem,  the  recurrence  relations  (3.26a,  b,  c)  hold  among 


the  coefficients  D  ”  (i?  ),  with 
r,  t  s 


fV 


./i?  cos0  i  1 
Jf  s  o 


(i?  icos0  ) 


(t-f): 


tV  i£  _cos0 
\  s  o 


A“fi(iC  )  = 

t  3 


(2t  +  1)  —  - 

(t+f):  Al-tVp/M+t  3\  t.  .  ’ 


M + 1  even 


M  + 1  odd 


(5.2) 


M  > 

For  the  Neumann  problem  the  coefficients  D  ‘  (i£  )  are  related  through  equations 

t  S 

(4.  22a,  b,  c)  with 
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M  ,  {iK tcosd  )M_1  n, 


MP 


.  /  i£  cos  9  t  r 
it  s  o 


t  V  i?  t  cos  9 

N  8  O 


.  P  /i?  cos0  rr 

—  tCi€  cose  V 

!  ,  s  o  t\  1?  ZcosO 
+  1  •—  x  s  o 


/i?  cos0  ±1\1 

-<t+*Xi?  tcosQ  )V  -j---  S  - 
S  O  t-lV  1?  tcosd 


M  +  t  even. 


o 

(5.3) 


and 


AtM,i(i?  )  =  0  . 

I  8 


M+t  odd 


(5.4) 


The  prime  on  Q  (i?  )  in  (5.3)  implies  differentiation  with  respect  to  if  . 

1  S  £  s 

The  Legendre  function  Q  (i?,)  in  (5.1)  must  now  be  redefined  since  can 

r  l  1 

now  assume  values  between  0  and  1  as  well  as  values  greater  than  1 .  This  has 
been  done  in  Appendix  H  where  we  show  that 

1 


'?  rjm 

Qm(  =  (-l)m  r(n  +  m+l)ni/2)  (r+1) 

"  '  in+12m  ffa+f)  (j+  /?Prm+1 


^F^n-m+l.  |-m;  n+|;  - 


1 


2  ’  9  2 

<$+£+!)> 


.  5  £0 


(5.5) 


The  far  field  is  given  by  (5. 1)  by  letting  5  — >  op. 
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uSf(Pj) 


ikcc,  a)  n  ...  .n-M 

1“ 4'("kC)  (n-M): 

1  n-0  M-0 


M 


J ..  ..XI.  I 


'  / _ j  /_ _ I>*  ’  (i?  )P  (n,)cosi0 

t=0  1=0  °’  s  1  1  1 


(5.  G) 


The  scattering  cross  section  is  given  by  (3.30)  and,  in  the  present  case,  can 
be  written  in  the  form 


OP 

i  2  *>n  \2n 
a  -  4tc  /  (kc) 


2n 


n=0 


m=0 


/  , ,  n-rm  sf  .  .  sf.  . 

(-1;  u  (P,)u  (p  )  , 

2n-m  1  m  1 


(5.7) 


where  we  have  taken  k  to  be  real,  and 

n  _  n-M  _ t 

tn-io:  • 


M-0 


Nose-on  Incidence 


When  0  =  0,  we  can  rewrite  (5.1)  as  follows: 
o 


(5.8) 


ike?. 


oo 


n  ...  _  .  M  M 

\  ’  +  V  K — "'v 

u  (Pi)  .  e  Zj-O' ’  -fe2^4-^,<«.WPi<V  • 


M-0 


M 


rVT 


(5.9) 


where,  for  the  Dirichlet  problem. 
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M+l..,.  .  2 

r,t  UV  ~  rlr+l)-t(t  +  l) 


“2  2 

r  _M  ...  .  _  t  ...  . 

2r  -  1  Dr-l,t  1?s)  +  2t  —  1  Dr,t-l(l?s) 

+  (t+1)2  M  ...  .  (r+1)2  M  .  .  r/t 
'  2t  +  3  Dr,t+i  1?s  '  2r  +  3  Dr+l,t1?sJ  ’  M  =  0,1,2,... 


(5.10a) 


A  ,  ...  . 

„M+1, ._  .  V  Qr  V  „M+1,..  4  .  .  M+l, ..  . . 

Dt|t  “«SI=-Z_,  W7Dr,ttl5S)+At  C5S);M  =0.1,2,.. 

r=0  t  s 


.  (5.10b) 


D°  (i?  )  =  A°(i?  )  , 
o,  o  s  os 


(5.10c) 


AX> 


(2t+  lHif  i  1)M 
s 


,  M  +  t  even 


0  ,  M  +  t  odd 


(5.11) 


The  corresponding  expressions  for  the  Neumann  problem  are, 


M+l, ._  . _ 2 _  r2  M  ...  .  _  t2  M  ...  . 

°r,t  1?s  “  r(r+l)-t(t  +  l)  |_2r- 1  °r-l,t(l?s)  +  2t  -  1  °r,  t-1  1?s 

+  (t+1)2  M  ....  (r  +  1)2  M  ...  r^t 
'  2t  +  3  r,  t+1  s'  2r  +  3  r+l,t  s  ’  M  =  0,1,2,... 


(5.12a) 


M+l...  , 
Dt,t  (l?s)  = 


V'W  M+l ...  uy  WnM  ...  .  .M+l, ._ . 

4*  Qt(iCs)’Dr,t(lV  r,t(lV  At  (1V' 


Q  (15.) 


M  =  0,1,2,...  (5.12b) 


an  =  o , 

o,  o  s 


(5.12c) 
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AX> 


-\[w  M 


t2t+l)(i?  ±1) 
8 


,  M  + 1  even 


0  ,  M  +  t  odd 


(5.13) 


In  both  the  Dlrichlet  and  Neumann  problems  the  far  field  is  given  by 


ikcf,  oo 


sf  e 

U  (P,»  =  — 


>Pt<v . 


1  n-0 


(5.14) 


and  the  scattering  cross  section  be  equation  (5.  7)  where  in  the  present  ca  e, 


(5.15) 


5.2  The  Disc 

As  we  mentioned  earlier,  when  ?  =  0  the  oblate  spheroid  degenerates  to  a 

s 

disc  of  radius  c  in  the  x,  y  plane,  with  center  at  the  origin  (Morse  and  Feshbach, 
1953,  p.  1292).  It  is  easy  to  verify  from  the  corresponding  formulas  for  the  oblate 
spheroid  that  the  scattered  field  due  to  the  presence  of  the  disc  is  given  by 


...  co  n  .  .n-M  M  M  t 

"V  *  e  XJ  ->ni>  X2^Z]Dr!'t,<C) 


•Q^i^lP^n^cosi , 


(5.16) 
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where,  for  the  Dirichlet  problem,  the  recurrence  relations  (3.26a,  b,  c)  hold  among 
M  i 

the  coefficients  D  '  (0),  with 
r,  t 


r 


Af'(0,  = 


3 


■f 


(tcos0  ) 


M 


1 


.M+l 


(2t+ 1) 


(t  -  i ) : 
(t 


o  , 


M  +  t  odd 


tV.cos  0  / 
o 


+,,:  Wy- 


(0) 


M  +  t  even 


(5.17) 


For  the  Neumann  problem,  the  coefficients  D  ’  (0)  are  related  through  equations 

r,  t 

(4.22a,  b,  c)  with 


A^'*(0)  =  cos0 


(t  cos  9  ) 
o 


M-l 


o  2M+1 


(t-f): 


W)'-  W*  !)«>>■  U*IK 


(M-t)(t-<  +  1)P 


t+1 


\pos0  J 


+  (M  +  t  +  l)(t  -f)P 


f  P-N 

t-l\COS0  / 


O  -J 


M  +  t  even  , 


AtM*  i(0)  =  0  , 


M  +  t  odd 


From  (5.16)  the  far  field  is 


ikcSj  _oo 

•  *7 


1  n=0  M=0  V  t=0  1=0 


n-M  M  t 


(5.18) 


(5. 19) 


(0) 


P  (q  )cosf0  . 


(5.20) 
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The  scattering  cross  section  is  given  by  (5.  7)  with 


n  ^  ^  ^n-M  M  t 

“fv  =  (-«"•'  -(„.■«)■  YU2<-»t>'  d”;'(0)p!(0)cos/)( 

M-0  '  t=0  1-0  0.1.1  1 


(5. 21) 


Normal  Incidence : 

When  6^  =  0  ,  we  can  write,  as  we  did  in  section  5.1, 


s,  .  lkc5l 
u  <Pj )  -  e 


_ |  n  ^  ^  ^n-M  M  ^  M 

jL<-lkc,n  X}-»“  -^mf  5S  D”  «(0)Qr(15l»Pt( 

n-0  M=0  r=0  t=0 


(5.22) 


where,  for  the  Dirichlet  problem,  the  coefficients  D  (0)  are  related  through 

r,t 

equations  (5. 10a,  b,  c),  with 


AtM(0)  = 


./ — , +  „.M  _ 2t+ 1 _ 

“  =»*■(»), r(JB.|)s,„ 


,  M+t  even 


0  ,  M  +  t  odd 


(5.23) 


IV1 

For  the  Neumann  problem  the  coefficients  D  X0)  in  (5.23)  are  related  thrc.igh 

r,  t  ^ 

(5.12a,  b,  c),  with 


a“(0)  = 


''  -V7(t 


1)M_I  M 


,  M+t  even 


0  ,  M  +  t  odd  . 
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VI 

NUMERICAL  CALCULATIONS 

As  a  demonstration  of  their  usefulness,  the  theoretical  results  have  been 
employed  to  calculate  the  scattering  cross  section  of  some  representative  prolate 
spheroids  for  both  Dirichlet  and  Neumann  boundary  conditions.  The  prolate  spher¬ 
oids  considered  had  major  to  minor  axis  ratios  of  10:1,  5:1,  and  2:1.  Back  scat¬ 
tered  and  forward  scattered  cross  sections  were  determined  as  functions  of  wave¬ 
length,  and  complete  polar  diagrams  of  bistatic  cross  section  were  obtained  for  a 
few  special  values  of  kc.  All  calculations  were  carried  out  for  a  plane  wave  incident 
along  the  axis  of  symmetry  of  the  spheroid. 

The  expressions  employed  for  this  calculation,  which  we  repeat  here  for  con¬ 
venience,  were  equation  (3.32) 


o  =  4rc 


and  equation  (3.38) 


oo  2n 

£  u®)2”  il  <- 


n=0 


m=0 


.n+m  sf  .  .  sf, 

11  u2n-m(Vum(V 


(6.1) 


.  .m-j 

^s4-^  i 

— ~r~  V  >P.(n,) 

(m-j)'.  o,  l  s  l  1 


(6.2) 


where  D*  .  in  (6.2)  is  given  by  (3.35)  for  the  Dirichlet  problem  and  (4.25)  for  the 

o,  1 

Neumann  problem. 

The  series  in  (6. 1)  was  terminated  at  n  =  10  for  the  2:1  and  5:1  spheroids 

and  at  n  =  9  for  the  10:1  spheroid.  Thus  the  cross  section  values  included  terms 

20  1 8 

up  to  and  including  (kc)  and  (kc)  respectively.  The  back  and  forward  scattering 
results  were  also  obtained  for  smaller  values  of  n  so  as  to  reveal  the  manner  in 
which  the  inclusion  of  higher  order  terms  improves  the  Rayleigh  approximation. 
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Figures  6-1  and  6-2  present  the  back  scattering  (r^  =  1)  cross  sections  of  soft 
and  hard  spheroids  respectively.  The  cross  section  values  are  normalized  with  re¬ 
spect  to  the  geometric  optics  value 


c 

g-O. 


=  nc 


2  2 
2  <5s‘  11 


(6.3) 


The  number  associated  with  each  curve  indicates  the  value  of  n  at  which  equation 

(6. 1)  was  terminated.  The  Rayleigh  curve  (the  curve  obtained  by  terminating  (6. 1) 

at  the  first  nonvanishing  power  of  kc)  is  denoted  by  n  =  0  for  the  soft  spheroid  and 

by  n  =  2  for  the  hard  spheroid.  The  exact  result  shown  in  Figures  6-1  and  6-2  was 

computed  from  the  prolate  spheroidal  function  series  (Senior,  1966).  Also  included 

in  each  figure  is  the  maximum  value  of  ka  (  =  kc?  )  for  which  the  series  in  (6. 1)  con- 

s 

verges,  i.e.  the  radius  of  convergence,  as  estimated  by  Darling  and  Senior  (1965). 
The  present  low'  frequency  calculations  have  no  precedent  except  in  the  case  of  the 
10:1  hard  spheroid  where  similar  calculations  (though  not  as  extensive)  were  re¬ 
ported  by  Sleator  (1964). 

Figures  6-3  and  6-4  present  the  forward  scattering  (r^  =-l)  cross  sections  oi 
the  same  spheroids.  The  cross  section  values  are  normalized  with  respect  to  the 
limiting  form  of  the  bistatic  geometric  optics  value 


2  2C2 

a  =  7ra  =  nc  ”  (6. 4) 

g.o.  s 

As  before,  the  number  associated  with  each  curve  designates  the  value  of  n  at  which 
the  series  in  equation  (6. 1)  was  terminated.  No  exact  results  were  available  for 
comparison  in  this  case. 

Figures  6-5  through  6-9  present  polar  diagrams  of  the  bistatic  cross  sections 
of  the  same  spheroids.  Since  the  polar  diagram  is  symmetric  for  nose-on  incidence. 


45 


THE  UNIVERSITY  OF  MICHIGAN 

7133-5-T 


which  is  the  only  case  considered,  each  figure  includes  data  for  both  hard  and  soft 
spheroids.  The  back  and  forward  scattered  cross  sections  lie  on  the  h^avy  vertical 
line  bisecting  the  figure  with  the  back  scattering  (0^  =0)  value  on  the  upper  part  and 
the  forward  scattering  (0^  =j)  value  on  the  lower.  The  values  of  the  cross  section 
are  normalized  with  respect  to  the  geometric  optics  cross  section,  viz.. 


a 

go. 


,  .4  2 
-  4jrb  a 


f' 


a2(l  +  n, 


)  +  b2(l 


v] 


-2 


=  4jtc 


2 

<*r1+V2 


(6.5) 


with  ijj  =  cos  0  ^  . 

As  noted  previously,  the  values  presented  for  the  2:1  and  5:1  spheroids  were  ob¬ 
tained  after  terminating  the  series  in  (6. 1)  at  n  =  10  while  for  the  10:1  spheroid  the 
series  was  terminated  at  n  =  9. 

Similar  calculations  have  been  carried  out  by  Spence  and  Granger  (1951)  for 

hard  spheroids  though  the  values  of  ?  and  kc  were  different  from  those  employed 

s 

here.  In  the  few  cases  where  comparison  was  possible  (kc  =  l,  a/b  =  5, 10),  good 
agreement  was  obtained. 


I 
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FIG.  6-2b:  BACK  SCATTERING  CROSS  SECTION  OF  HARD,  5:1 
PROLATE  SPHEROID  FOR  NOSE-ON  INCIDENCE. 


51 


THE  UNIVERSITY  OF  MICHIGAN 

7133-5-T 


FIG.  6-2c:  BACK  SCATTERING  CROSS  SECTION  OF  HARD,  10:1 
PROLATE  SPHEROID  FOR  NOSE-ON  INCIDENCE. 
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FIG.  6 -3 a:  FORWARD  SCATTERING  CROSS  SECTION  OF  SOFT,  2:1 
PROLATE  SPHEROID  FOR  NOSE-ON  INCIDENCE. 
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FIG.  6-3b:  FORWARD  SCATTERING  CROSS  SECTION  OF  SOFT,  5:1 
PROLATE  SPHEROID  FOR  NOSE-ON  INCIDENCE. 
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FIG.  6-3c:  FORWARD  SCATTERING  CjIOSS  SECTION  OF  SOFT.  10:1 
PROLATE  SPHEROID  FOR  iIOSE-ON  INCIDENCE. 
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FIG.  6.4a:  FORWARD  SCATTERING  CROSS  SECTION  OF  HARD,  2:1 
PROLATE  SPHEROID  FOR  NOSE-ON  INCIDENCE. 

56 


2 

7r  a 


-8 

10 


0  .2  .4  .G  .8  1.0  1.2  1.4  1.6  1.8  2 

ka 

FIG.  6- 4b:  FORWARD  SCATTERING  CROSS  SECTION  OF  HARD,  5:1 
PROLATE  SPHEROID  FOR  NOSE-ON  INCIDENCE. 
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FIG.  6-4c:  FORWARD  SCATTERING  CROSS  SECTION  OF  HARD,  10:1 
PROLATE  SPHEROID  FOR  NOSE-ON  INCIDENCE. 
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FIG.  6-5:  BISTATIC  CROSS  SECTION  OF  2:1  PROLATE  SPHEROID 
FOR  NOSE-ON  INCIDENCE. 
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FIG.  (>-(>:  BISTATIC  CROSS  SECTION  OF  5:1  PROLATE  SPHEROID  FOR 
NOSE-ON  INCIDENCE. 
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FIG.  6-7:  BISTATIC  CROSS  SECTION  OF  5:1  PROLATE  SPHEROID 
FOR  NOSE-ON  INCIDENCE. 
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FIG.  6-8:  BISTATIC  CROSS  SECTION  OF  10:1  PROLATE  SPHEROID 
FOR  NOSE-ON  INCIDENCE. 
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FIG.  6-9:  BISTATIC  CROSS  SECTION  OF  10:1  PROLATE  SPHEROID 
FOR  NOSE-ON  INCIDENCE. 
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APPENDIX  A  + 

THE  REGULARITY  OF  THE  FUNCTION  u(p)  =  e“lkC  ?  “  n  u(p) 

In  this  appendix  we  offer  a  proof  that  the  function  u(p)  in  (2. 14)  is  regular 
at  infinity  in  the  sense  of  Kellogg,  that  is 


hni  j  rUp)  j  <  on  and 
r  — >  oo 


lim 
r  — *  ao 


2  aup) 

r  3r 


<  oo  , 


0^:9  7r 
0^^2tt  . 

(A.  1) 


The  proof  is  based  on  an  expansion  theorem  (Wilcox,  1956)  which  guarantees  ihat 
ihe  field  scai.ered  by  the  prolate  spheroid  may  be  written  in  the  form 

ikrx“ o  f  (0,0) 

uS(p)  =  —  ,  r  >  a  (A.  2) 

n=0  r 


where  the  series  is  uniformly  and  absolutely  convergent  for  all  r,  0,0  provided 
r  >a,  a  being  the  radius  of  the  smallest  sphere  completely  enclosing  the  prolate 
spheroid. 

s 

From  (A.  2)  it  is  clear  that  u  (p)  satisfies  the  first  of  conditions  (A.  1)  but 
not  the  second  and,  consequently, is  not  regular  at  infinity.  The  function  u(p),  how¬ 
ever,  which  by  (2.14)  and  (A.  2)  may  be  written 


oo 


uXp)  =  e 


-ik(cc;  -  r icq)  1 


f  (0,0) 
n 


r  ‘ n 
n=0  r 


(A. 3) 


can  be  shown  to  satisfy  the  Kellogg  conditions.  The  proof  is  as  follows: 

The  variables  £  and  rj  are  related  to  the  spheroidal  coordinates  by  the 

T 


equations  1 

?  "  2c 


FT  2  [2  S 

y r  +2crcos0  +  c  +  y r  -2crcos0  +  c 


■iff 


2  2 
+  2crcos0+c 


-{F 


2cr  cos0  +  c 


n  2c  LV 

The  factor  c(?iq)  appearing  in  the  exponential  of  (  A. 3)  can  now  be  written 
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Employing  (A.  5)  and  (A.  9)  in  (A.  7)  we  have  that 


lim  r 


2  3cj(p)  I 


r->oo 


e*  ikccosfl  +  OO/r)  ^ 

|_  n=0  r11  * 


+  ikcos0. 

=  e  f  <  oo  , 

o 


(A.  10) 


which  shows  that  the  second  Kellogg  condition  holds  also. 
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APPENDIX  B 

THE  SURFACE  INTEGRAL  FOR  THE  DLRICHLET  PROBLEM 

In  this  appendix  vve  evaluate  the  surface  integral  of  (3.4).  Repeating  the  ex¬ 
pression, 

IS(Pj)  =  -c(-^  -  1)  j  ^  p  )  .  (B.l) 

From  equations  (3.2)  and  (3.3), 


-ike(-  _ /])  -ikc(cosO  ?  i j  +  sind  l/£  -1  41-))  cos  0) 

,  .  s  os  o1  s  '  . „ 

Up  )  =  -e  e  (B. 2) 


and  from  (2.21),  with  £  <£ 

s  1 


00  n 


-  G  (p. .  p  )  T — /  ,/  ,(-D  f  (2n+l)  7— — -  cosm<0-0) 

d£  o  1  s  — -j  r)  m  (n+m)'.  r  rl 

s  nr0  m=0  L 


.-|2 


Pm(o,)Pm(o) 

n  1  n 


P‘V  ) 

p"V  )'Qm(g)---  S-  Qm(g  )'QnV.) 

n  s  n  1  nni(f  )  n  s  n  1 
m  '  s7 


oo  n 


4ttc(?j  -  1)  n=0  m  0 
s 


c  (2n  +  1)  f~~rr  cos  t.i(0 -  0  ) 
J  ^  m  (n  +  m).  r  rl 


QmU.) 

p"V  >?"’<„>  -2-*- 

n  1  n  ,m 


(B.  3) 


where,  above,  we  used  the  Wronskian  relation 
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[~Qm(a  pm(?) 

L  n  n  _ 


Q  (?  )  ’ 

n  s 


,m 


=  Pm(?)'Qm(?)-Qm(€)'Pm(?)  =  —p—  • 

n  n  n  n  1 


IB. 4) 
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Substituting  (B.2)  and  (B.  3)  in  the  integral  (B.  1),  we  have 


-ike?  ct>  n 

s  e  - 

I  <p,>  =  -Z-7- 

*  n-0  m-0 


y  >  .  ,  v(n-m):  ^m,  ,  n  1 

/  ,/  c  (2n+l)  — — -  P  (nJ - 

r1^ — ^  m  (n  +  m)'.  n  1  ^m.. 


W 


-ikc(?  cos 6  tl)r| 
\  .  so 

*\  dne 


s2  77 


P™<n)  \  d0e 

Jo 


i-i 


r 2  2 

-ike sin0Q  y'?g- 1  yl  -  r)  eos0 


cos  m(0-  0^). 
(B.  5) 


The  functions  involved  in  the  integrands  are  continuous  in  the  intervals  of  integration 
and  the  only  assumption  we  made  in  interchanging  integration  and  summation  is  the 
uniformity  of  convergence  of  the  series  (cf.  Whittaker  and  Watson,  1952,  p.  78). 

We  now  use  the  expansion  (Magnus  and  Oberhettinger,  1949,  p.  155) 


oo 

ikp  cos  0  \  .m  _  ,,  ,  , 

e  =  /  l  e  J  tkp)cosm0  . 

z ^  mm 
m=0 

r 2  /  2 

Utilizing  (B. 6)  in  (B. 5),  with  kp  =  -kcsin@  \?-l  yl-n  ,  results 

o  s  ' 
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-ikc(?  cos0  tl)n 
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-IKCVC  COS  U  -  Iff]  r~z -  - 

drje  "  °  Pm(r])Jf(kcsin0oy?s- 1  yl J 


M 


r)  y\  d0cosm(0- 0^)cosf0  . 

0 

(B.7) 


The  integration  with  respect  to  0  can  be  simply  performed,  while  to  integrate  with 
respect  to  rj  we  use  the  relation  (Morse  and  Feshbach,  1953,  p.  1325), 
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,  lz  cos  vcosn  _m.  »  , 

d^e  P  (cosp)J  (z  sini.'sinu)sinu 

n  m 


.n-m  j  2  z 


P  (cosidJ 
n  n+ 


In  this  expression  we  let 

n  =  cos p 

zcost'  =  -kc(F  cos0  tl) 
s  o 

f~2 

z  sini'  =  kcsinO  >/?  -  1 
o  V  s 


so  that 


=  [kc(r  -  COS  0  )~| 
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S  O 
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We  can  then  write 
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jVvIk 
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s  o  J 


Performing  the  0  integration  in  (B.7)  and  using  (B.9),  we  obtain 
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where,  above,  we  used  the  relation  (Magnus  and  Oberhettinger,  1949,  p.  (3) 


m  .n+m  _m.  . 

P  (-x)  =  (-1)  P  (xj  » 
n  n 


-1  ^  x  $  1 


(B.  1 1) 


We  now  expand  the  Bessel  function  in  (B.  10)  according  to  (Magnus  and  Ober¬ 
hettinger,  1949,  p.  16), 


(iz/2)2* 
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(B.  12) 


to  get 
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-ike? 
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2  M^O  S  °  m  (f+m)! 


n  cose  ii\ 
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(B.  13) 


where  in  the  above  series  in  k  the  only  nonzero  contributions  to  the  coefficients  are 
made  by  terms  for  which  M  +  f  is  an  even  integer. 

We  have  then  written  the  surface  integral  (B.  1)  as  a  power  series  in  kc  of 
the  form 


IS(p)  =  e  ?s  r.i-iko^) 

1  fco  M  1 


(B.  14) 


where  I^ip^)  is  given  by 
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APPENDIX  C 


,cx> 


EVALUATION  OF  THE  INTEGRAL  \  C  (?,?  ?  )Q  (?)d?  . 

\  t  l  s  r 


'00 


d?Ct<?,?l,?s)Qr(?)  =  \  d?Qt<?1)Pt(?)Q!-(?)  +  \  d?Pt(?l)Qt(?,Qir(?) 
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From  Legendre's  associated  equation  we  have  that 


_d 

d? 


_d 

d? 


-  ,  dQ '(?)' 


,  <!</(?) 

<1-5  )-4— 

d? 


f2  1  i 

rlr+J)--—  !  Q  (?)  -  0 
■i  i  r 


1  -? 


t(t+l)  - 


1-rJ  1 


Q  (?)  =0 


t  i 

Multiplying  the  first  of  these  equations  by  Q  (?)  and  the  t.eeond  by  Q  (?)  and  sub- 

t  r 

traeting  the  second  from  the  first  we  obtain  the  following: 
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(1-0 
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d? 
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Integrating  this  expression  we  have  that 
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Similarly, 
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Hr+l)-t(t  +  l)  d?1  ‘W  d?1 


rir+i)-«t+l)  Q<(5s)  iw  dis  dis 


*,V  f>, 
"tt:  1Qt< 


Qr(V  (t  +  w  l  W  .  (t+0:  ,  J 

r<r+ 1) -t(t+ 1)  (t -i)l  '  '  r(r+l)-t(t  +  l)  l._  ,  (t-i)l  '  ; 

VV 
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APPENDIX  D 

DERIVATION  OF  RELATION  (3.21) 

In  this  appendix  we  give  a  proof  of  the  statement  of  equation  (3.21),  that  is 


>M;V. . 

r,  t  s 


DTfS  )  =  (±1) 


r+t_M,  f._  . 

Dt  (5)  ,  f.r.t^M. 

i,  r  s 


(D.  1) 


The  choice  of  sign  is  determined  by  the  sign  chosen  in  (2. 18). 

The  proof  follows  an  inductive  argument.  First  we  show  that  (D.  1)  is  true 
for  M  =  0  and  then  that  if  it  is  true  for  any  M  it  is  true  for  M  +  I. 

Since  r  =  t  =  0  when  M  =  0,  equation  (D.  1)  is  certainly  true  for  M  =  0. 
Assume  next  that  it  is  true  for  M.  We  can  then  integrate  (3. 19)  and,  following  the 
same  procedure  as  we  did  there,  end  up  with  the  recurrence  relations  (3.26a,  b,  c). 
We  are  interested  mainly  in  (3.26a)  since  for  r  =  t  equation  (D.  1)  is  obviously  true 
Repeating  here  (3.26a)  and  subsequently  employing  it  in  (D.l)  which  is  assumed  to 
hold  for  M,  we  obtain 


dm+1,<(s  )  = 


r,  t  s  r(r+l)-t(t+l) 


rtr-O  M.l  .  Kt-f)  M,< 

2r-l  Dr-l,tV+  2t-l  Dr,t-l(5S’ 


+  (t+l)(t  +  i  +  l)  M,  f  (r+l)(r  +  f  + 1)  M,  t 

2t+  3  r.t+l'V"  2r+3  r+l,tV 


i^1f-(tl)I'ft"1DM’i(?) 

2r-l  t,r-l  1 


r(r  +  l)  -  t(t+  1) 


t(t-l)  +  rft-lnM  J  .  +  (t  +  l)(tj J±l),  +  ur+i+lMt 
+  2t-  1  “  1  Dt-1,  r  s  ”  2t  +  3  (~}  Dt+l,r?s 


(r+l)(r+C  +  i)  ■  r+t+1  M.l 

2r +  3  ’  t.ru'V 


(cont'd) 
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=  (tl) 


t(t+l)-r(r+l) 

+  (rjJ.Mr+f  +  1)  M,  i 


r^>D- 1 

[_2t-  1  t-1,  r 


(€  >*  ) 

s  2r- 1  t,  r-1  s 


2r+3 


Ws>- 


(t-HMt+i  +  1) 
2t  +  3 


d**;*  (?  >] 

t+l,r  s  | 


=  (Il) 


t,r 


V 


So  (D.  1)  is  true  for  M+ 1  if  it  is  true  for  M.  Since  it  is  true  for  M  = 
for  all  M  (M  =  0,1,2, ...). 


(D.  2) 

it  is  true 
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APPENDIX  E 

DERIVATION  OF  EQUATIONS  (3. 24a,  b) 

In  order  to  arrive  at  equations  (3.24a,  b),  we  start  with  (3.  22)  and  (3.23) 
which  we  repeat  here 


M+l  M+l 


rM+l(p,) 


i 


t-0  r=0  f=0 


r(r+l)-t(t  +  l) 


M,l 

2r- 1  r-l,tV 


.  t(t-/)  M,/  ,  +  (t+l)(t  +  f +  1)  ^M, t  ..  . 

4  Dr,  t-l(V  -  ~ 2TT5 - Dr,  t+l(V 


(r  +  l)(r  +  /  +  l)  M,  i 

2r+3  r+l.t'V 


Q*<€J  -  . 

“^■Qt(€l)"Qr(5l) 


JLW 


P^rj^cosf^  , 
(E.  1) 


M+l  M+l  .  t 


'm+M  VQr<5l)P!(Vcos,f!l  .  (E.2) 

t=0  r=0  1=0  '  brill  l 


When  r^t,  a  comparison  of  these  two  equations  gives  (3.24a).  When  r  =  t,  we  re 
write  the  above  equations  as  follows. 

Equation  (E.  1): 


M+l  M+l  .  t 


IM+ltpl)  2 


1 


t=0  r=0  1=0 


r(r+  l)-t(t  +  l) 


(C ) 

2r-l  r-l.t'V 


-,L(t-OnM,f  (  +  (t+l)(t+i  +  l)  M,f 

+  2t  —  1  Dr,t-l(V-  2t  +  3  Dr,t+lUs) 


(r+lHr  +  f  +  1)  M,  l 
2r  +  3  r+1, t  s* 


Qr^I^Pt^nl)COsf  ^1 


M+l  M+l  .  t 


-2 


- - J  - 1  r(r  +  +  1) 

t=0  r=0  1= 0  1  ;  u  ; 


ITT  Dr  \  ,«„> 

2r-  1  r-1,  t  s 


(cont'd) 
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-  t(t-i)  nM ,t  tm  + 
+  ITT  Dr,  t-lls8>  ' 


(r+  l)(r+f  + 1) 
2r+3 


(t+l)(t+f +  1)  M,f  .  . 
2t+  3  r,t+lVV 


-j  ) 


W 


Equation  (E.2): 


M+l  t 


kw  ■  2  2 


A  comparison  of  (E.  3)  and  (E.  4)  gives 


M+l 

t’t  8  ^ 


(E.3) 


(E.4) 


-2 

r(r+l)-t(t+l) 


[_2r-l  Dr-l.tV+  2t-l  Dr,t-l'V 


+  (t  +  l)(t+i  +  l)  M, l  .  .  (r+lKr+l  +  1) 

2t+3  rft+lVV”  2r+3 


Using  (3.24a),  the  above  expression  can  be  written 


M+l.  ,  f.  . 

8  TTTn  /->z /e  \  r»l  8 


^M,f  71  Qr 

j, 


<4<V 


w 


(E.5) 


(E.6) 


r=0  Qt(?s) 


which  is  equation  (3.24b). 
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APPENDIX  F 

THE  SURFACE  INTEGRAL  FOR  THE  NEUMANN  PROBLEM 
The  integral  to  be  evaluated  in  this  appendix  is  that  of  equation  (4 . 6) 

+1  r2n 


s  2  'ikC?« 

Kp,)=-c<?  -l)e  s 

]  s 


.  _  .  >  :  3U  (P  ) 

dr)  \  d0e+ikcr,G  (p  ,  p) 

\  o  1  s 


3? 


(F.l) 


'-1  JO 


By  equation  (3.3) 


3u*(p  ) 
s 


=  -ike 


cost?  n+sin0  r 
o  o  / 


2 

cos0 


v* 


u(ps) 


r 


=  -ikcjcos0  rj+siriO  ~~r  n  1 

l  »  »7?r 


^s/l  “02  ~1  -ikc|cos  0o?sr)  +  sin0Q^f^-  ljl  -  r)2  cos^J 

■  cos  0  L 


1 


Substituting  the  above  expression  together  with  the  appropriate  part  (S  <  £^)  of 
(2.22)  in  (F.  1)  we  have 


!%,)  = 


-^(?2-l)e  ?SX^Z](~ire  (2n+1) 

47 r  s  — rri  in 

n=0  m=0 


,  m 


(n-  in)'. 


Pm<?  >’  -  r+1 

Pm(?  )  -  ~  ~  Qm(€  ) 

_n  c  Qm,?  }I  n  s 

ns  0-1 


dr)  Pin(rj)  e 
n 


(n  +  m)'. 

-ikcl?  cos0  t]"]n 
is  o  -> 


n  l  n  l 


'2*  r  r  1  2 

d0  cos  6 ^n  + sin 9 ^ -]===== — cos  f 


)0 


i  2  r 2  2 

y  1  -  n  ~ |  -ike  sin0Q  y?g  -  1  yl-r)  cosf) 


°/ 


9 

S  -1 

S 


X  cos  m(0-  0^)  .  (F.  2) 
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£  Gin0  'r — i  \ — i  .  ,  q  ip  ) 

ike  s  o  \  \  .  ..m+1,  .  (n-m)'.  m,  ,  n  V 

M.-i  W 


cd  n 


m.  . 

(n-m)'.  m,  4  VV 


-ikc(£  cos  6  t  l)n  . - *  x  —  - - -> 

dne  0  yl  -n  Ucsmojt*- 1  /l- n2 


(F.4) 


To  perform  the  integration  with  respect  to  n  we  employ  the  following  recurrence 
relations  (Magnus  and  Oberhettinger,  1949,  p,  62) 

(2n+l)r)Pn  (n)  =  (n-m  +  l)P^(n)  +  (n  +  m)P™  ^(r)),  m  £  n,  n  =  0, 1,  ...  (F.5) 

(2n+l)  P™(r?)  =  P™  ^(ri)~  P^^n).  m^n,  n  =  0,1,2,...  (F.6) 

(2n+l)  Jl  -  rj  P  (r?)  =  (n- m+ l)(n  -  m +  2)P™  /(r))-(n-m- l)(n  +  m)Pm  /(n), 
*  n  n+1  n-1 

m  ^  n,  n  =  0, 1,  . . .  (F.7) 

Substitution  of  these  expressions  in  (F.4)  and  a  simple  rearrangement  of  the  terms 


leads  to 


ikc£ 

s  rs.  .  ike  . 
e  I  (p  )  =  —  cos0 

Y  c  O 


oo  n 


;,m  J~(n-  m)’.  nm  .  ,  ^n-i^l^ 


n=0  m=0 


/  J  114/.  nUI  ,  ^ 

e  (-1)  <  7— — rr  P  ,(rj,)- 

»  m  (n  +  m)l  n-1  1 


Qm  .(£  )’ 

n-1  s 


+  lco.n.0, 

<0+m)-  ”+1  1  <}"«■  )J  1 

n+1  s 


-ikc(£  cos0  tl)n 


s  o  „m 


P>Wm(csi“9„^  a’1 


(cont'd) 
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(n-m+1)1.  m  ,  .  (n-m  - 1)1  m  .  Si-1^1^ 

(n  +  m+1)'.  n+1 ’V  m  ^  (n+m- 1)1  n-i  V  m  ^ 

“  V*n+i  ?s;  ^n-l  *s' 


r  oo  n+1 

.  ike  n  S  \  \ 

+  -r-  sinO  /  / 

0  1  n=6J  m=0 


<-«"  ./^  j  1((Z)pm-‘( 

V  z  n+V2  n 


^  <p~  m+ 1)1  m  .  .  Si-1^1^  (n-  m-H)'.  m  .  .  ^n+1^1^ 
I  (n+m- 1)'.  n-r'V  m  .  'n+m- 1)1  n+lV  m  . 

L  VlV  Qn+l(V 


cos 


where 


z  =  kc(£  tcos0  ) 
s  o 


cos  m0j , 


(F.8) 


(F.9) 


?  COS0  tl 
S  0 

?  +CCS0 
S  0 


(F.10 


Equation  (F.8)  is  now  put  in  the  following  form 

.  r-  oo  n 

ike  Jyv  .  ..n+1 


m  <C<V 

P  (n,) - cosm^ 

"  1  Qm(f  )•  1 

n  s 


a>  n 


l/2^",  ,  »~m  (n-m+1)! 

VT  V3/2(Z)PZ+1(|3)  (n+m)'. 


Qm(C,) 


m  >z  n  ^  "-1  (n+m-1>1  n  1  Qm(5  )■ 

n  s 


#  cosn^ 


(cont'd) 
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p  r 00  p~2 

m 

•P  (n,)-— - cos  m0 

n  1  Qm(C  >’  1 

a  s 


J  1((z)Pm^(f9^SlI 

y  z  n-y2  n-1  M  (n+m)'. 


oo  a 


n=U  m=C 


J'AjL  t  (n- m):  m.  . 

I  *  V%(z>Pn+1  W  (oTrnK  P„  <V  “STT 


r  r-  00  n+2 

,  ike  .  „  \  '  \  ' 

H - —  sin0  •  /  > 

2  O./XT  rJ 

VV 1 


in  ■  hj  a  i 

n+l  ^  (n+m)l  VV^^,  cos 

n 


|<-i j*lJ&  ljKV1?"*!211 

1  v  z  n+-y2  n+l  H  (n+m)'. 


p“<  >  Q"m<5i  * 

“  ‘  Wv  c°sm0‘ 

n  s 


2S“-‘ffv.4»c: 


(n-m)l  m,  v  Qn  V 

(p>  H+rr,  9\.  p  'V - 

(n  +  m-  2)'.  n  1  m( 

n  s 


cos  mi 


Substitution  of  the  relations  (Magnus  and  Oberhettinger,  1949,  p.  16) 

'-1! ;«-TrV%w+£v%« 


(F.  11) 


(F.  12) 


Tn+3/2 


(F.  13) 


in  (F.  11)  and  a  regrouping  of  the  terms  leads  to 
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ikc5s,s,  ,  t  ,E 
e  I(Pl)  =  -Tco«  flj- 


go  n 


,n=0  m*0 


,  ,4n  2n+l  (n-m)'.  „  ,  , 

e  (_1)  ~ o —  7~~ — TT  J  ,  i/(z) 
m  2z  (n  +  m)'.  n+72 


n  Qm(?J 

•  (n+m)Pm  (^-(n-m  +  DP11!  .(,3)  PRn,) — ” —  cos  m0 
n-l  K  n+1  J  n  1  m(  rl 


Q  15  )' 
n  s 


ao  □ 


R  Re  (-i)n(2n+l)(3Pra(/3)J  „(z)' ro-^.  P°V)--°  *■-  cosm(l 

m  n  n+%  (n+  “>•  "  1  qP'is  )■  1 

*n  3  o  1 


00  n 


-if  Sine  -pZ-fe 
2 


.n=l  m 


,  ,  ..c(2n+l)  ,  v  (n  m)l 

1;  2z  n+V2  /J  (n+m)l 


_  Qm(F  ) 

Pm^1(j3)  +  P^t^iS)  Pm(n, )  — — ~  cos  m 0 
[_n-l  11+1  °  1  Qm/c  v  1 


Q  (?_)’ 

□  s 


00  n 


(n-m)'. 


2mm  OV 

z  „m.  „m,  .  n  1 


4  Wi<«'  ^ 

n=l  m=l  ^ 


Qm(?  )' 
n  s 


00  n 


^  \  \  /  ,vn  (2n+ 1) 

cosm^- 2-J-')  ~2T~ 

n=l  m  =  l 


,  4  Rn  -  m  +  2)'.  „m-l,„4 
V%(Z)L  <n+mK  -+1  ® 


,  (n-  m) 

(n+  rn  -  ! 


nr.  _m-l,  l4  „m,  4  1  a 

— 7^r,  P  ,  ((3)  P  <r?, )  - cos  mp 

-2)1  n-l  !J  n  1  m(f  ,,  '1 

ns 


(F.  14) 


This  expression  can  be  simplified  using  the  properties  of  the  Legendre  functions 

mentioned  above.  After  simplifying  and  collecting  terms  in  J  x/(z)  and  its  deriv- 

n+  I2 

ative,  we  have 
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ikC?sTs,  .  kc  [2tt 

e  1  (pi}  =  "T  v~ 


r-  cd  n 


.  n=0  m=0 


e  (~i)P(2n+l)  J  .!  ,(z)' 

m  (n+m)!  n+Y2 


5  i/l -(3“  Qm<S.) 

p cosy  +  -f  ■  —  sinO  Pm(,i)Pm(n,)  — ~ - cos  m0  ■ 

L  °  /?rr  si  “  "  1  Qm(s  r  ‘j 


W 


i-  CD  n 


kc  /  2ir 
2  V  z 


.n=u  m=c 


,  ,  ,vn  2n  + 1  (n-m):  „  ,  v  „  ,  , 

,e  (-1)  — ^ —  7— 7 — cose  (n+m) P  .(j 3) 
— ^  m  2z  (n  +  m)'.  n+Y2  o  n-1 


-cose  (n-m+DP™  ((3) -sine  -7  f  —  |/l  -|}2  Pm(0) 
o  n+1  o  /  2  *  n 

rV1 


+  sine  /3Pm  1(^)-(n+m)(n- m+l)sine  pr-  -  ]3  Pm  *(jJ) 

-  °W'  *  - 


—Dii  .  n  1  j 

P  (r>  ) -  cosm0. 

“  QX>' 


(F.  15 


But  by  (F.  10) 


?slM 

8  cose  +  i  =■ - -  sine  =  1 

0  te7'  ° 


(F.  16 


Moreover,  from  the  definition  of  /3  and  the  recurrence  relations  (Magnus  and  Ober- 
hettinger,  1949,  p.  62) 


Pm+1(|3)  =  -p===  (n-m^P^Mn  +  nOP™.^)  , 
Vl-0  L  -1 


(F.17; 
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Pm-\p)  = 
n 


1 


E 


(n+m)  fl  -]:2  L  “ 


P>-Ci(<3) 


we  have  that 


cos 0  (n  +  m)Pm  (p)-cos0  (n- m  +  l)Pm  (3)  -  sin.6  ___ 
o  n-1  o  n+1  o  f_ 2 


Pm(^) 

n 


+  sin0 


“pJ=^pm  *(p) -(n  +  m)(n- m  + l)sin0  <3 Pm  ^ 

Xi  n  °P7lp  n 

1  S  |  s 


03) 


±(?  i  cos  0  ) 
S 


-  tV3  V  /  p  _ ^ 

~2 - “  (n-m+l)p“  (0)*  ?  P™(3)-(n+m)Pm,(.W  I 

2  L  n  *  1  s  n  fl“i  J 


Substituting  (F.  16)  and  (F.19)  in  (F.  15)  we  have 
ike? 


oo  n 


2U  z_ 

'  nr0  m=0 


( 1 1)(?  ±  cos  0  ) 


oo  n 

kc  I2ir  '  ,  ,vn  £a+J.  (n-  m)'.  .  s  o 

2z  (n+m):  n+Vz  , 


(n-m  +  DP1"  (p)  +  ?  Pm(p)  -  (n  +  m)Pm  (p) 
n+1  s  n  n-1 


Qm(?,) 

„m.  .  n  1  u 

P  (rv) -  cos  me, 

"  1  Qm(f  )'  1 

n  s 


88 


(F.  18) 


(F.19) 


(F.  20) 
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A  simple  inspection  cf  this  expression  reveals  that  it  is  zero  for  n-:0.  We 
can  therefore  write 

co 

..  .M  s  .  . 

(F.21) 


-ike? 

.s,  „  s 


Hpj)  =  e  “  ^_J(-ikc)Ml^I(p1) 
M=i 


where 


ts  in  \  ic  +, 

2J 


W =  VC°8<V 


M  < 

M-l  NT"  V" 

/ 

m 


2f+ 1 


(i- m)l 


(T)(?  I  cos  9  )  r 

(2M+l)pJn(/3)  + - T - 2 

C-i 

a 


* m  (¥>  ,<+m,! 
(f-m+i)P^1(3)  +  ?sPjn(^) 


Qm(?  ) 

-(i+mjP^1  (/3)  I  ^  P!"(n-) — ~ — ~  cosmfl  ,  M  +  f  even  , 
f-!  J  /  i  m.,  1 


1  1  q  “<?_>• 

b 


(F.  22a) 


I-^p.)  =  0  ,  M  +  i  odd  , 

M  1 


(F.  22b) 


where  above  we  have  substituted  (F.9)  for  z  and  we  have  rearranged  the  series. 
Equation  (F.22a)  can  be  further  simplified  by  taking  into  consideration  (F.  10)  for 
P  and  the  relation  (F.  5).  In  this  way  we  can  write 


M  f 

Cv  ■  2  l>r<gp(m<v<<v-n*i 


(F.  23) 


1=0  m=0 


where 
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APPENDIX  G 

f°°  i  l 

THE  INTEGRAL  \  d?K  {?,?  ,?  )Q  (?) 

J?  1  1  8  r 

s 


According  to  the  results  of  Appendix  C, 


GOO  pCj  pCD 

^  dc  p'fSKi'ffl  +  Pp"^  dfQpKjp) 


W'  i 


'CO 


i  'V,l 

w 


0.(5,)  \  d5Q:(5)Q  (5) 


J? 


^t^l^  2  ft  l  £  I  ~ I 

fPuTTT)  V«1  -  1>lft«1)Qr€€1>* 


r(r+l) 


>  + 
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.A 
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APPENDIX  H 
REDEFINITION  OF  Q 


m 


m, 


The  original  definition  of  Q  (/i)  as  given  by  (2.26)  is 


,  (-l)m  r(n  +  m+l)r(l/2)  (i? -  l)m^2 

V"’  ‘  7o+i - r.  J.  1ST  2f 

2  /  (0+2)  M 


/n  +  m+_2  n  +  m  + 1  ^  lA 

IV  2  ’  2  ;n  2;  2 )’ 


|/i|>l.  | arg(/i -  1) |  <  Jr . 


(H.l) 


Hobson  (1953,  pp  233-234)  has  shown  that  if 


z  =  *1+  hi  -  1  , 


(H.  2) 


then  the  function 


1 


,  2  2 

u(M)  =  -  1)-- 

M  n+m+1  2 

z 


Fi(!+n,,n  +  m-H;n+f;i)  , 


I Z  |  >  1 ,  |arg(M  -  I)")  <  7r 


(H.  3) 


satisfies  the  associated  Legendre  equation.  Using  this  expression  we  can  define  a 


new  function  Q m(iu)  which  holds  for  jzj  >  1  cr  equivalently  |^|  >  0,  which  is  iden¬ 


tical  to  Q™(/i)  givenj?y  (H.l)  in  their  common  domain  of  definition,  |/i  j  >  1 .  To  do 
this  it  is  sufficient  to  compare  (H.  1)  and  (H.3)  for  ’arge  values  of  \n\ .  The  result¬ 
ing  relation  between  the  two  functions  is 


QmWM-l>ro2m-(a+m+‘>r(l/a  „(„) 

"  r<n+f> 


(H.4) 


or 


94 


THE  UNIVERSITY  OF  MICHIGAN 

7133-5-T 


2  2 m 

Qm(  )  =  (  Dm 2m_P(n+m+l)T(l/2)  (^  -1)  F 
n  f(n+|)  zn+m+1  2 


I(i+m.D+m+l;n+!ii), 


|z|>l,  |arg(ji- l)j  <  n.  (H.5) 


Letting  ^  =  i?,  ?  >  0,  we  have  z  =  i  (?+  y?  + 1  J ,  and 

,_2  .m/2 

nm  (-2)m  r(n+m-4)r(l/2)  ' 

QnU5)  '  “5T  r(n+|) 


f  i  3 

x  2FlV 2  +  m'  n  *‘m  +  1;n+  2;  ■ 


•  ffnj- 


?  >0  . 


(H.6) 


Using  the  relation  (Magnus  and  Oberhettinger,  1949,  p.  8) 


2F1(a,  b;c;z)  =(l-z)C  a  ^F^c-a.c-bjcjz) 


we  can  write 


_2m  s 

-F1(i+m'“  +  m+l;n+|;-4)  -  2F1("-m+1'i-m'n+!^) 


(H.7) 


Letting  z  =  iU+j£  +1  )  and  substituting  in  (H.6)  we  obtain 


(-2)  r(n+m  +  l)r(l/2) 

Qn  (1?)=“E+T  r.  3. 

i  f(n+-) 


_ (g2+l)m/2 _ 

_(?+  l/?^)2  +  l]  (?+  1)°' 


,-3m+l 


'2Fl(n_m  +  1'2'm;n  +  !; 


l/?n)5 


,  no  (H.  8) 


7 133- 5- T 


which  holds  at  ?  =  0  also  since  (n-m+l)  +  (~-m)-(n  +  ^)  -  -2m  £  0  for 
ir  =0,1,2,...  (Magnus  and  Obcrhettinger,  1949,  p.  7). 

Equation  (H.  8)  can  be  rewritten  to  read 


(~l)m  r(n+m  +  1 )  r  ( 1  /  2) 


Qm(i?)  = 
n 


(?2+d 


-  “  m 


i"+l  2m 


r(n+2) 


12  An-m+1 

+o 

V 


2Fi(n-m  +  1,|-m;n  +  |; 


or; 
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APPENDIX  I 

THE  FAR  FIELD  FOR  THE  DISC 

In  this  appendix  we  give  the  first  six  terms  in  the  far  field  expansion  for  the 
disc  for  both  Dirichlet  and  Neumann  boundary  conditions  with  the  incident  wave  at 
normal  incidence. 

The  far  field  in  both  cases  is  given  by  (5.23).  This  expansion  involves  the 
M 

coefficients  A^  (0)  given  by  (5.23)  and  (5.24)  which  contain  the  Legendre  functions 

of  the  second  kind  and  their  first  derivatives  with  respect  to  i£  evaluated  at  £  =0. 

s  s 

Their  values  are  determined  as  follows. 

From  equation  (5.5), 


Q  lii:)  =  -L  r(n+nr(i/2)  1 

"  r1  r<n+f>  (c^/12'1 


Letting  ?  =  0,  we  have 


F1(n+1,-;n+~;. 


?  £0 


(I.D 


1  r(n+l)r(l/2)  „  ,  .  .  1  .  3 

Qn(0)  -  ^Tl  - “ sT"  2Fl(n+1'  n+ 2;'1} 

1  r(n  +  ~) 


(1.2) 


Now, 


2Fi(n+l,  2‘  n+  2’ 


J7_ 

„n+l 


<n  +  | 


r(f+1)r(!+1 


(1.3) 


Then  (1.2)  becomes 


Q  (0)  = 

n 


jrnl 


(1.4) 


(2i 


n+1 


Lr(l  +  1)J 


:,s  See,  for  example,  Handbook  of  Mathematical  Functions ,  National  Bureau  of 
Standards,  Applied  Math.  Seri  T'o.  55,  p.  557  (June  19o4). 
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Or,  we  can  write 


Q2n0)  ,„..2n+l  ’  n 


Q2n+1(0)  = 


(2i)  n'.n! 

ff(2n  +  i)'. 


(-4) 


n+i  r„,  ^ 

lr{n+2U 


,  n  =  0, 1, . . . 


(1.5) 


(1.6) 


Turning  now  to  the  derivative  of  Q  ,  denote  by  Q'(0)  the  derivative  of  Q  (if)  with 

n  n  n 

respect  to  if  evaluated  at  f  =  0.  From  (1. 1) 

VU)  ■  7  ~9~;1)r<3l/2>  {-<n+l>2F1(nM. n+f;-l) 

l  r  vn+  — )  *- 


2(n  +  1) 
2n+  3  2 


F1(D  +  2,  2;n+2;_1^} 


(1.7) 


Employing  (1.3),  (1.4)  and  the  relatk 


F^a,  b;  a- b; -1)  =  2  a  fn  (b-  1)  1  P(a-b  +  2) 


-1 


2  1 


1 


Ln|a,r(f +|a-b) 


1 


n2  +  2a)r(1+2a'b)- 


(1.8) 


which  can  be  found  in  the  same  reference  and  page  as  (1.3),  equation  (1.7)  becomes 

n+1  2 


Q'(0)  =  i(n+l)Q  (0)  -  T 
n  n  ,nnn+l 

1  ^ 


Ln|  +  i)r,-"  +  „  r(s+I)r(|+|)J 


.  (1.9) 


Finally,  with  the  help  of  (1.5)  and  (1.6),  we  obtain  for  (1.9), 
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Q2n(0)  = 


Jr  (2n)! 


c 

(-4)nfr(n+^)] 

~  _J 


n  =  0, 1, 


(1. 10) 


(2i)  n!n'. 


(1.11) 


We  now  employ  equation  (5.26)  and  write: 
For  the  Dirichlet  Case 


u8f(p  )  =  -  P  (rj  ) 

0  1  7T  O  1 


“iV  -  ^  W 

71 


U2(V  9tt  P2(r,l)+V  3  ~97r)Po(V 

X7T  ' 


Sf. 

U3  <Pl>  = 


Sf.  V 

U4  lpl> 


sf. 

W 


w+(gy,oi)w+(3j- g3\ 

iW+(-5+i)p2(v 


+ wJ  P„(V 


~^~2  P4ini)+  V  ~4  +  ~~2)  P2(  V 
525^  ^  9*  567*  '  1  1 


f 64  80  ,  508  >\  _  ,  , 

\~6  -  ~+ - 2j  P„(V 

\  9i  2025)7  y 


Substituting  these  results  in  (5.25)  we  obtain 
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Sf,  . 
u  <PL> 


7~  i'f  p„<V+ikc'lp<.S,+k2c2  -£w 

1  l__  r  •— 

h(4 - £) W]  +  ik3°3  [_a  p2(r>i,+  (- -1  +  72) p„<  V 

JT  ->  1—  V  7T  '  —I 

^  M5x  P4<V+(^+T^jW+(‘ai  +  73'^)Po<,’l) 


7T  tiff 


+  ^  ^2P4(V+(-71+^-2)  P2(V 


Y64  80  .  508  \  „  ,  . 

+(  -7  - — 7+ - -  P  (n. ) 

9r  20257:  ^  ° 


P  <n.)  +0(k6c6)l  . 

y  °  1  >  . 


For  the  Neumann  Case 


uSf(p  )  =  0 
o  1 


u*f(pi)  =  0 


U2f(pl)  =  -'^P1(’>1) 


uf<Pi)  -  0 


U  In  )  —  -  P  (n  )+ -  P  (n  ) 

4  -1  75tt  3  V  75tt  l'V 


(1.12) 


sf.  .  4  „  .  . 

u=lpi)  =  -riW 

2m 
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